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Abstract 

Host-parasite co-evolution can have profound impacts on a wide range of ecological and evolutionary pro¬ 
cesses, including population dynamics, the maintenance of genetic diversity, and the evolution of recombi¬ 
nation. To examine the co-evolution of quantitative traits in hosts and parasites, we present and study a 
co-evolutionary model of a social parasite-host system that incorporates (1) ecological dynamics that feed 
back into their co-evolutionary outcomes; (2) variation in whether the parasite is obligate or facultative; and 
(3) Holling Type II functional responses between host and parasite, which are particularly suitable for social 
parasites that face time costs for host location and its social manipulation. We perform local and global 
analyses for the co-evolutionary model and the corresponding ecological model. In the absence of evolution, 
our ecological model analyses imply that an extremely small value of the death rate for facultative social 
parasites, primarily due to hunting/searching for potential host species, can drive a host extinct globally 
under certain conditions, while an extremely large value of the death rate can drive the parasite extinct 
globally. The facultative parasite system can have one, two, or three interior equilibria, while the obligate 
parasite system can have either one or three interior equilibria. Multiple interior equilibria result in rich 
dynamics with multiple attractors. The ecological system, in particular, can exhibit bi-stability between 
the facultative-parasite-only equilibrium and the interior coexistence equilibrium when it has two interior 
equilibria. Our analysis on the co-evolutionary model provides important insights on how co-evolution can 
change the ecological and evolutionary outcomes of host-parasite interactions. Our findings suggest that: 

(a) The host and parasite can select different strategies that may result in local extinction of one species. 
These strategies can have convergence stability (CS), but may not be evolutionary stable strategies (ESS); 

(b) The host and its facultative (or obligate) parasite can have ESS that drive the host (or the obligate 
parasite) extinct locally; (c) Trait functions play an important role in the CS of both boundary and interior 
equilibria, as well as their ESS; and (d) A small variance in the trait difference that measures parasitism 
efficiency can destabilize the co-evolutionary system, and generate evolutionary arms-race dynamics with 
different host-parasite fluctuating patterns. 

Keywords: Bistability, Evolutionary Game Theory, Evolutionary Stable Strategy, Holling Type II 
Functional Response, Facultative/obligate Parasite. 


1. Introduction 

Parasitism takes many forms in nature, but can be broadly defined as a symbiosis in which one member 
(the parasite) benefits from the use of resources gathered by the other member (the host). In most cases, the 
parasite lives on and/or feeds directly from the host. In social parasitism, however, the parasite manipulates 
its host behaviorally rather than physiologically, and derives benefit from work provided by the host (Cini 
et al 2015). The different forms of social parasitism exemplify the diversity of ways in which social parasites 


1 Sciences and Mathematics Faculty, College of Letters and Sciences, Arizona State University, Mesa, AZ 85212, USA 

(yun.kangOasu.edu) 

2 School of Life Sciences, Arizona State University, Tempe, AZ 85281, USA (j .fewell@asu.edu) 


Preprint submitted to XXX 


May 19, 2015 





can cheat by manipulating the efforts of their hosts. In brood parasitism, a parent deceives or dominates 
other individuals into rearing their young (Field 1992; Brandt et al 2005; Spottiswoode et al 2012). For 
avian brood parasites, this generally occurs via deception, when females lay eggs on other birds’s nests and 
the host parents respond as if the parasitic chick is one of their own (reviewed by Kruger 2007; Davies 2011). 
This deception is often continued by the parasitic chick, which acts as one of the host brood but may claim 
a disproportionate proportion of the food resources provisioned by the host parent (Kilner et al 2004), or 
even remove the hosts own chicks from the nest (Spottiswoode et al 2012). 

In another taxonomic realm, insect social parasites, or inquilines, gain acceptance into a social insect 
colony by mimicking the pheromonal signals of queens or workers; They then live within that colony, ex¬ 
ploiting its resources for their own reproduction (Bourke and Franks 1991; Buschinger 2009). Colony para¬ 
sitism can also occur via more direct aggression, as when queens of some wasp species directly displace other 
queens and usurp the colony as their own (Field 1992; Shorter and Tibbets 2009). In brood raiding ants, 
colonies invade nearby nests and retrieve brood, usually in pupal form just before they morph into new adult 
workers. The emerging workers behave as offspring of the invading colony, performing the same functions 
as if they were in their natal nest (Brandt et al 2005; Holldobler and Wilson 1990 &2009). In some of the 
most extreme examples of brood raiding, ant species (slave-making ants) become dependent on the raided 
workers, without which they cannot functionally maintain their own colonies. Their raiding efficiency and 
intensity has co-evolutionary consequences for both the host and raiding species (Foitzik et al 2001). 

Although diverse in their specifics, each of these cases takes a temporal progression in which the parasite 
must locate a host and then manipulate it into social acceptance of either the parasite or its offspring; it does 
so either by deception or aggression. It then establishes an ongoing relationship in which the host provides 
care, usually to the parasite’s offspring, by providing resources, defending the parasitic offspring, and/or 
providing other parental care. These stages provide a scaffold around which we can model the interplay of 
ecological and evolutionary effects on parasites and their host. The dynamics within the social parasite-host 
system can also provide a framework to capture the dynamics inherent in parasite-host relationships more 
generally. 

Interspecific social parasitism, in which one species parasitizes another, is rarer than intraspecific social 
parasitism, but can have profound effects on community dynamics. Its impact on the relationship between 
host and parasite is highly dependent on host number (Sorenson 1997; Brandt et al 2005). When modeling 
the dyadic relationship between a parasitic species and a given host, we must consider whether that parasite 
is fully dependent on that specific host, i.e. obligate, or whether it can utilize other host species as well. A 
parasite with multiple potential hosts (generalist parasites), essentially behaves as if it is facultative within 
the context of that dyadic relationship. For the purposes of this model, we focus on interspecific parasitism, 
and facultative and generalist parasites are considered as equivalent. 

The ecological and evolutionary dynamics for facultative versus obligate social parasites are very differ¬ 
ent. A generalist parasite, such as the brown headed cowbird may parasitize tens to hundreds of species, and 
must flexibly deal with the associated variation in parental care and diets across those species (Rothstein 
1975). In contrast, the common cuckoo specializes on a few species, focusing primarily on the reed warbler. 
The reed warbler has evolved defensive strategies including recognition of cuckoos and their eggs (Davies 
and Brook 1988; Stoddard and Stevens 2011). In return, the cuckoo uses a variety of mimicking strategies, 
from chick calls to egg color and markings to overcome the reed warbler’s defenses (Davies and Brook 1988; 
Kilner et al 1991; Davies 2011; Stoddard and Stevens 2011). The ecological and evolutionary drivers of 
obligate versus facultative parasitism are complex. Phylogenetic analyses suggest the number of hosts for 
cowbirds has increasingly expanded evolutionarily (Rothstein et al 2002); however, this parasite is also in¬ 
creasing in range and number. In contrast, the cuckoo host-parasite relationship suggests no clear pattern 
of expansion over evolutionary time, and variation in host number may be better explained by ecological 
conditions (Rothstein et al 2002). 
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Obligate social parasitism generates a potentially tight dynamical fitness relationship between parasite 
and host. As the host suffers fitness costs from the parasite, to avoid extinction it must evolve defensive 
strategies to counter the parasite (Bogusch et al 2006). In turn, the parasite, dependent on the resources 
acquired from its host, is selected to overcome the host’s defensive strategies (Poulin et al. 2000). As a result, 
the continuous interactions between the parasite and its host lead to co-evolutionary dynamics, as illustrated 
by the evolutionary arms-race paradigm that has been used for many host-parasite systems (Anderson and 
May 1982; Thompson and Burdon 1992; Foitzik et al 2003). Close co-evolutionary interactions in a stepwise 
fashion are especially likely to occur when host and parasite exhibit similar generation times and population 
sizes; this is generally the case for social parasites, because brood parasitism requires a match with host 
offspring developmental timing (Foitzik et al 2003). In such situations, hosts are expected to more closely 
match the parasite in strategy evolution, and thus to evolve increased resistance strategies when parasite 
pressure is strong (Foitzik et al 2001& 2003). In support, recent studies on co-evolution in both avian and 
insect social parasites found indications of arms races and resistance strategies in highly parasitized host 
populations (Foitzik et al. 2001; Kilner and Langmore 2011). 

The evolutionary dynamics of host-parasite interactions can be influenced by multiple intersecting pa¬ 
rameters. These may include: the efficiency of parasitism on a specific host; the genetic structure of host and 
parasite populations; migration rates of parasite and host, and; the degree of mutual specialization (Brandt 
et al 2005). The dynamics among these variables can be complex, and their intensity may vary considerably 
between systems for which the parasite is obligatorily bound to a single host, versus those for which the 
parasite is facultative on a given host and/or can exploit multiple hosts. In this paper, we develop a simple 
co-evolutionary model by using evolutionary game theory (EGT), to investigate the ways in which these 
co-evolutionary dynamics can change the ecological and evolutionary outcomes of host-parasite symbioses. 

Surprisingly, the mathematical models that examine the co-evolution of quantitative traits in hosts and 
parasites are relatively few in number (Hocherg and van Baalen 1998; van Baalen 1998; Gandon et al. 2002; 
Koella and Boete 2003; Restif and Koella 2003; Bonds 2006; Zu et al 2007; Best et al 2009). However, 
they provide an important theoretical framework within which to consider the importance of co-evolutionary 
dynamics to the evolution of hosts and parasites. Hocliberg and van Baalen (1998) used simulations to 
study host-parasite co-evolution in relation to spatial heterogeneities. Restif and Koella (2003) showed that 
when host defense occurs specifically through avoidance, the host-parasite relationship can move towards a 
single co-evolutionary stable state (CoESS), which cannot be invaded by other strategies. Van Baalen (1998) 
showed that resistance through recovery, rather than avoidance, can generate bistability such that resistance 
and virulence are either both low or both high. 

A limitation of the above models is that they examined only the evolutionary stability (ES) of the 
outcomes, i.e. their ability to be invaded once reached. Eshel (1983) and Christiansen (1991) examined con¬ 
vergence stability (CS), the other important component of the evolutionary process, determining whether 
populations will move away from a CoESS over evolutionary time and whether evolutionary branching will 
occur. Dieckmann and Law (1996), and Marrow et al. (1996) provided a general understanding of the impact 
of co-evolution on both ES and CS. In another approach, the work of Kisdi (2006) studied the impact of 
trade- offs on co-evolutionary dynamics through a general model of predator-prey co-evolution. Zu et al 
(2007) expanded this to examine evolutionary dynamics of the host (as prey) and the parasite (as predator) 
with a Holling Type II function. Their results showed that branching in the host can induce secondary 
branching in the parasite, and that these evolutionary dynamics can produce a stable limit cycle. The work 
of Best et al. (2009) also emphasized the importance of considering co-evolutionary dynamics and showed 
that certain highly virulent parasites may result from responses to host evolution. 

These models collectively assume that the parasite (or predator) is obligate. However, facultative and/or 
generalist social parasites are actually the norm. Further, many models assume a slow evolutionary timescale 
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to allow application of a timescale separation argument and reduce dimensions and simplify mathematical 
analysis (Feng et al 2004). However, in social parasitism the changes in allele frequencies (and associated 
phenotypes) likely occur at the same rate as changes in population densities or spatial distributions, which 
can alter the ecological processes driving changes in population densities or distribution (Gingerich 2009; 
Cortez and Weitz 2014). Motivated by these issues, we present a fully co-evolutionary model of a social 
host-parasite system that includes consideration of: (1) the ecological dynamics that feed back into the 
co-evolutionary outcome; and (2) variation in the form of social parasitism from obligate to facultative. The 
interaction between host and parasite is modeled by using Holling Type II functional responses. These are 
particularly suitable for brood parasitism, because they allow consideration of the need for parasites to search 
for and locate hosts, as well as time interacting with hosts. Our model allow us to explore the following 
questions: 

1. How do ecological dynamics change when a parasite transitions between being facultative and obligate? 

2. Can the parasite or host have evolutionary stable strategies (ESS) that can drive each other extinct at 
local or global scales? 

3. Can co-evolution rescue a host-parasite symbiosis from extinction? 

4. What are the effects of trait functions on the ecological and evolutionary dynamics between the parasite 
and host? 

The remaining sections of this article are organized as follows: In Section 2, we provide the background 
of the EGT modeling approach, define CS and ESS mathematically, and derive our co-evolutionary host- 
parasite model to incorporate a parasite can be either obligate or facultative. In Section 3, we perform and 
compare local and global analyses of the ecological dynamics of the co-evolutionary model in the absence 
of evolution. Our results show that the ecological model of parasite being facultative can generate one, two 
or three interior equilibria with the possibility that the host goes to extinction either locally or globally. 
However, the obligate parasite model can have either one or three interior equilibria with the possibility 
that the obligate parasite goes extinct. In Section 4, we study the dynamics and related ESS of the fully 
co-evolutionary model of the host-parasite system, including ecological dynamics that feed back into the co- 
evolutionary outcomes. In particular, when trait functions follow normal distributions, we derive sufficient 
conditions that the host-only equilibrium and the parasite-only equilibrium can have ESS, and study the 
boundary and interior dynamics. The work of the latter case shows that: (a) the co-evolutionary model 
can have multiple interior equilibria where each equilibrium potentially has two ESS; (b) evolution can save 
the host from extinction; and (c) a small variance in parasitism efficiency can destabilize the system, thus 
generating evolutionary arms-race dynamics with different host-parasite fluctuating patterns. In Section 5, 
we conclude our study by providing a summary of our results and their potential biological implications. 
The last section gives the detailed proofs of our analytical results. 


2. Evolutionary game theory models: co-evolutionary host-parasite models 

In an EGT model, the co-evolution of hosts and their parasites can be considered as a mathematical 
game, where the host and parasite are players, with corresponding strategies, strategy sets, and pay-offs 
(Vincent and Brown 2005). The strategies represent phenotypic traits with heritable components, and the 
strategy set is the collection of all evolutionarily feasible strategies for a particular individual, with corre¬ 
sponding fitness pay-offs. Here we use EGT to investigate host-parasite co-evolution within a population 
dynamical framework, under the assumptions that the interaction between a host and parasite is modeled 
by a Holling Type II functional response, and that the reliance of the parasite on the identified host can 
range from completely obligate (the parasite is completely dependent on that host for any fitness pay-off) 
to facultative (the parasite can survive without parasitism and/or can parasitize other available hosts). 
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2.1. The modeling framework 

We follow the modeling methodology for EGT presented in Vincent and Brown (2005), (also see references 
Abrams et al 1993a, 1993b; Rael et al 2011; Cushing and Hudson 2012; Kang and Udiani 2014; Kang et al 
2015). The methodology derives equations that describe the population dynamics of n interacting species 
Xi together with the dynamics of (mean) phenotypic traits Ui (or strategies) which serve to characterize 
all individuals of a particular species and are assumed to have a heritable component. Let x = [si,£ 2 ] 
denote the vector of population densities of host, parasite and u = [mi, 1x2] denote the vector of all strategies 
u = [iii,M 2 ] G [Ui,U 2 ] = U used by the species i, which are distinct and drawn from the same set Ui of 
evolutionarily feasible traits. The host-parasite co-evolution dynamic models in terms of differential equations 
have the following form: 


dt — Xi(~*i {Vi , 3 ?) | y i= Ui — XiHi(u,x) 

duj _ 2 dGj(vj,u,x) I (1) 

dt i dvi \vi=Ui' ) 


where G,(m,,m, x) is the fitness of a focal individual that chooses (or inherits) trait v,- L when the population 
has mean trait u and density x\ Hi(u,x) is considered as the fitness function for species Xf. and of is the 
variance in traits (strategies) present in species Xi about the mean trait ut which therefore a measure of the 
“evolutionary speed”. The host-parasite population dynamics and its associated trait dynamics (1) together 
constitute a dynamical system for the 4-vector 

[x,u] = [xi, £ 2 , mi, M 2 ] € R+ x Ui x U 2 

whose dynamic describes an evolutionary process known as Darwinian dynamics (Vincent and Brown 2005). 
These equations (1) allow us to investigate the role of evolution plays in determining, for example, the dy¬ 
namical outcomes of host-parasite interactions. 

Darwinian dynamics will often possess a (locally asymptotically) stable equilibrium [x*,m*]. Suppose 
this equilibrium remains stable when additional species (with their associated traits) are added to the com¬ 
munity. This means that the coalition of traits Ui associated with those species present in the equilibrium 
(i.e. for which x* > 0) is able to resist invasion by other species with their other traits. For example, the 
4-dimensional equilibrium [x*,u*] of (1), when embedded in the higher 2m-dimensional space of the larger 
community (m > 2), remains (locally asymptotically) stable. In this case, we say that this coalition of traits 
u* is an evolutionary stable strategy (ESS). The ESS maximum principle provides a necessary condition for 
a coalition of traits in an equilibrium [x*,u*] to be an ESS (Vincent and Brown 2005). According to this 
principle, if a coalition of strategies is an ESS, then the G-function evaluated at equilibrium conditions, i.e. 
G(v,u* ,x*), takes on an isolated global maximum with respect to v at each trait Vi = u* in the coalition. 
Furthermore,this maximum value must equal 0 for the differential equation model. 

We say an equilibrium x* is an Ecologically Stable Equilibrium (ESE) if x* is locally asymptotically stable 
for the ecological dynamics ^ = x,Gj(Mj, m, x)\ Vi — Ui = XiHi(u,x) for given values of traits Mj. We say an 
equilibrium (x*,m*) has Convergency Stability if it is locally asymptotically stable in the co-evolutionary 
model (1). We say an equilibrium (x*,u*) has u* as an Evolutionary Stable Strategy (ESS) of the co¬ 
evolutionary model (1) if (x*,u*) has CS, and it satisfies the following ESS maximum principle: 

max{Gi(uj, m*, £*)} = G,(m*, u*, x*) = 0 for both * = 1,2 when x* > 0, i = 1,2 or (2) 

max{Gi(Mi, m*, x*)} = Gi(u*,u*,x*) < 0 for both * = 1,2 when x\xz = 0 (3) 

•u*EU* 

where u* = (m*, m^), x* = (x*, £ 3 ) an d Ui is the feasible set of trait values m which can be a bounded set or R. 
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2.2. Derivation of a host-parasite co-evolutionary model 

We consider a host-parasite co-evolutionary model with Holling Type II functional responses where the 
parasite can be facultative or obligate. Let Gi(vi,u, x) be the fitness of a focal individual that inherits trait 
Vi when the population has mean trait u and density x. We can take Gi in (1) with the following forms: 


Gi(vi,u 2 ,x) 
G 2 (v 2 , Ui,x) 


n(l 


Xi \ 
Ki(vx) ) 


a(v1 ,U2)X2 
l-\-h a(v i ,U2)xi 


e a(ui ,V2)xi 
l-\-h a(ui ,V 2 )xi 


d(v 2 ) +r 2 (l 


X 2 > 
-^ 2 (^ 2 ) ; 


( 4 ) 


where a(v 1 , 02 ) denotes the parasitism efficiency of a parasite with phenotypic trait v 2 on host individuals 
with phenotypic trait v\ with the assumption that the stronger host-parasite interactions are, the more 
similar host and parasite traits are; K-fvi) is the carrying capacity of species i individuals with phenotypic 
trait Vi] and d{v 2 ) > 0 is the death rate of parasitic individuals with phenotypic trait v 2 due to their hunt¬ 
ing or attacking all potential hosts. For simplicity, we assume that other parameters (i.e., the host intrinsic 
growth rate ry > 0, the parasite intrinsic growth rate r 2 > 0 in the absence of parasitism towards host species 
Xi, the parasite handling time h, and the conversion efficiency e) are not influenced by the quantitative traits. 


As a consequence, the co-evolutionary dynamics of mononrorphic resident host and parasite populations 
with traits u\ and u 2 is given by the following set of nonlinear equations: 


dx 1 
dt 


XiGi(vi,u,x)\ vl — ul = x 1 H 1 (u 1 ,u 2 ,x 1 ,x 2 ) 


Xl 




a(ui,U2)x2 
l-\-ha(ui ,U2)xi 


dx 2 
dt 


du\ 

dt 


x 2 G 2 (v 2 ,u,x)\ V2 = U2 = x 2 H 2 (u 1 ,u 2 ,x 1 ,x 2 ) 


X 2 


ea(u 1 ,u2)xx 
l-\-ha(ui ,U2)xi 


d + r2 ( x _ 


J dGi(vi,u,x ) | _ _2 dH 1 (u 1 ,u 2 ,x 1 ,x 2 ) _ 2 

1 dv i !«!=!!! 1 du\ 1 


rixi K[(ui) 
K i(ui )' 2 


*a' att( £' U2) 

,U2)xi)' 2 


( 5 ) 


dU 2 

dt 


^2 dG2(v2,U,x) | _ dH2(lli,U2,Xl,X2) _ ^2 

U Q Q - CLo Q - (Jr\ 

* OV 2 \V 2 = U 2 ^ OU 2 ^ 


9a(u 1 , 1 * 2 ) 
eXl du 2 _ 

(l+ha(ui,U2)xi)‘ 2 


d’{u 2 ) 


r 2 x 2 K!,\u2) 

K 2 (u 2 y 2 


where trait functions Ki(ui ), i = 1,2 and a(ui,u 2 ) are positive, bounded and smooth functions in Ui x U 2 . 
The co-evolutionary model (5) incorporates the following ecological assumptions: 


1. In the absence of other species, the population of each species follows a logistic growth function. When 
a host species x\ is absent, the population dynamics of the parasite x 2 depends on the death rate d 
caused by that parasite on all potential host species, coupled with the parasite’s intrinsic growth rate 
r 2 > 0, in the absence of parasitism on host species x±. For a chosen trait value u 2 of parasite x 2 , 
if d(u 2 ) > r 2 , then the parasite population with trait u 2 will go extinct in the absence of host ay. 
In such a case, we consider parasite x 2 to be obligate. If d{u 2 ) < r 2 , then the parasite population 
with trait u 2 can persist in the absence of host x±, and for purposes of this model is called facultative. 
Thus, in the case of obligate parasitism, the parasite species is unable to persist in the absence of the 
host. Under the case of facultative parasitism, the parasite can either (a) function in the absence of 
any parasitism (fitting the classic definition of facultatively parasitic), or they can survive using other 
host species. Assuming parasite-host dynamics are evolutionary labile, this formulation allows us to 
investigate conditions under which parasites might move evolutionarily from facultative to obligate 
and/or from specialist to generalist. 

2. The interactions between host and parasite can be described using a Holling Type II functional response 
equation to describe the average attacking rate of a parasite, time spent searching for and handling a 
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host, and time spent on other activities associated with the host-parasite relationship (Anderson and 
May 1978; Skalsii and Gilliam 2001).The Holling Type II functional response is particularly suitable 
for modeling social parasitism, because it allows consideration of the main challenges social parasites 
face: finding a host or host colony, overcoming the defenses of that host, and exploiting its resources for 
the parasite’s own reproduction (Cini et al 2011). As an example, in the case of avian brood parasites, 
overcoming defenses could include the tactics employed by a parasitic adult female and in relation to 
the host pair, while exploitation would consider the relationship between the chick and the host parents. 


3. In the absence of co-evolution, i.e., a\= a\ = 0, the co-evolutionary model (5) reduces to the following 
host-parasite ecological model: 


dx i 
dt 


Xi 




ax 2 

l-\-hax\ 


dx 2 
dt 




l-\-hax 




which can be host v.s. facultative or host v.s. obligate parasites depending on the the values of r 2 and 
d: parasite X 2 is facultative if r 2 < d while it is facultative if r 2 > d. 


4. Social parasitism also occurs facultatively within species. However, for purposes of this model, we 
assume that the host and parasite are different species; this allows us to consider them ecologically 
removed that there is negligible competition for food, space or other resources beyond the social par¬ 
asitism relationship. 


In the following section, we will explore the ecological dynamics of the co-evolutionary model (5) when 
of = of = 0 (in the absence of evolution). 


3. Ecological dynamics of a host-parasite model: facultative v.s. obligate parasites 

Let of = of = 0, then the co-evolutionary dynamics of model (5) is reduced to the following ecological 
model (6) with fixed traits Ui,i = 1,2 of host and parasite, respectively: 


= xiGx{v,u,x)\ v=ul = x 1 H 1 (u 1 ,u 2 ,x 1 ,x 2 ) = x 1 n (l - K ^ Ul ^ 


a(ui ,U2)X2 


l-\-ha(u\ ,U2)xi 

= X 2 G 2 {v,U,x )\ v = U2 = X2H 2 (u 1 ,U2,X 1 ,X 2 ) = X2 ~ ^ ( 1 ~ Ttffe 

whose state space is defined asT = {{xi,x 2 ) € R+}- In addition, we define the following two sets: 


( 6 ) 


dx 2 
dt 


X 1 = {(®!,0) G X} and A' 2 = {(0,x 2 ) G X}. 


For convenience, we let A'i(ui) = K\, K 2 (u 2 ) = K 2l a(u\,u 2 ) = a , and d(u 2 ) = d for short in this section. 
We have the following lemma regarding the basic dynamical properties of Model (6): 


Lemma 3.1 (Basic dynamical properties). Let K^ 1 = sup tligUi {A'i(ui)} and A 2 = sup U2gU2 {A' 2 (u 2 )}. 
Then the host-parasite model (6) is positively invariant and every trajectory starting in X attracts to the 
compact set C = [0, K^ 1 ] x [0, A ° 2( 'f —^j. 
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Notes: The results of Lemma 3.1 suggest that our co-evolutionary model (13) is biologically well-defined. 
The host-parasite model ( 6 ) always has the following two boundary equilibria E 00 = (0,0), E w = (A'i,0), 
and it has an additional boundary equilibrium Eq\ = ^0, (1 — f^)K 2 ^J on the y-axis if r 2 > d. The results 
on the boundary equilibria of Model ( 6 ) are summarized with the following proposition. 


Proposition 3.1 (Boundary equilibria). The host-parasite model ( 6 ) can have two or three boundary equi¬ 
libria where their existence and stability is listed in Table 1: 


Boundary Equilibria 

Stability Condition 

Eqo 

Saddle if r 2 < d; Source if r 2 > d. 

E'lo 

Saddle if r 2 + > d; Locally asymptotically stable if r 2 + b < d 

E'en 

Requires r 2 > d for the existence: it is a saddle if — > K 2 
while it is locally asymptotically stable if ^ < K 2 ^1 — T 

o 

A 

i—1 


Table 1: The local stability of boundary equilibria for Model (6) 


Notes: Proposition 3.1 implies that for the host-only equilibrium (A'i,0) to be locally stable (i.e. local 
extinction of parasite) requires that the death rate of the parasite due to searching for all potential host 
species must be larger than the sum of the parasite’s intrinsic growth rate r 2 in the absence of parasitism 
towards to Xi and benefits from the parasitism of x 1 , i.e. d > r 2 + ■ In this case, the parasite is 

obligate. However, for the parasite-only equilibrium ^ 0,/\2 ^1 — ^)) l° ca Hy stable (i.e. local extinc¬ 

tion of host) requires that the death rate of parasite due to searching for all potential host species must 
be smaller than the ratio of sum of the host’s intrinsic growth rate rq to the parasitism efficiency a, i.e. 
> K 2 ( 1 — ) >0. This is the case when the parasite is facultative. 


Theorem 3.1 (Persistence of the system). The host x\ of Model ( 6 ) is persistent in X if 


d > r 2 or — > K 2 ( 1 —— ) > 0. 


The parasite x 2 of Model (6) is persistent in X if 


eaK i 

r 2 + ——> d. 


1 + ahK\ 


Moreover, Model ( 6 ) is permanent in X if either 


eaKi 

r ‘ + T+mT, >d>r2 


or 



> 0 . 


Notes: Theorem 3.1 indicates that host x\ is always persistent if r 2 < d , while the parasite is always per¬ 
sistent if r 2 > d. The condition r 2 > d indicates that the parasite is facultative, i.e., it is able to survive in 
the absence of host x\. The coexistence of both host and parasite requires proper values of d regardless of 
the parasite being facultative or obligate. 












Let F(x i) = C 3 xf + c 2 xf + c\X\ + Cq where 


C 3 = rir 2 (ah) 2 

ci = a 2 hKiK- 2 (r 2 + f — d) + rrr 2 (l — 2ahK 1 ) 


c 2 = ahr x r 2 (2 — ahK\) 


c 0 = ar 2 K\ 



ri 

a 


and 

1 —C2 — \/c 2 ~ 3ciC 3 2 ~ c 2 + \/ c 2 ~ 3 ciC3 

*'c = - V - < 2 c = - V- • 

0C3 6 C 3 

Theorem 3.2 (Interior equilibria). Model (6) can have none, one, two or three interior equilibria de¬ 
pending on the values of parameters. The sufficient conditions for the existence of interior equilibria are 
listed in Table 2. Let (a:*, £ 2 ) an interior equilibrium of Model (6), then it is locally asymptotically stable 


Number of Interior 
Equilibria 

Sufficient Conditions 

None 

One of the following conditions holds: (1) d > r 2 + : or 

( 2 ) 7 T < K 2 (l — 7 ^) and 

a 2 r 2 K x K 2 (h + £ - > ^W+D 2 . or (3) K± < J_ and 

7T < K * (l - i) ; or (4) Kl > A. and r '£ a t“ h «f < 1 - 

One 

One of the following conditions holds: (l)r 2 < d < r 2 + ; 

or (2) ^ > k 2 (1 _ A) and (h + ± ; or (3) 

7T > Ki (l — 4), (h+ ± < r i a XlZ )2 and F(4) < 0. 

Two 

^<K 2 {l-±),{h + ±--^)< and F{xl) < 0 

with x 2 > 0. 

Three 

7T > K * (l - i) - { h + £ - 7?) < -^k\kT » F ( g c) > 0. and 
F(x 2 ) < 0 with x\ > 0. 


Table 2: The existence condition of interior equilibria for Model (6). 


if the following conditions are satisfied 


ri a rir 2 

— > — and - 

K 2 2 K x Ki 


+ 


(1 + ahK x y 


> a 2 hK 2 (r 2 + -). 


( 7 ) 


Notes: Theorem 3.2 indicates that the ecological model (6) can have one, two, or three interior equilibria 
when parasite is facultative, i.e., r 2 > d; while (6) can have one, or three interior equilibria when parasite 
is obligate, i.e., r 2 < d. Theoretically, it is possible that (6) has two interior equilibria in the singular cases 
when r 2 < d. We did not investigate such singular cases. The result for ignoring these singular cases is 
that host persists in Model (6) whenever r 2 < d. And this is different from the case when r 2 > d. The 
detailed results on the persistence have been provided in Theorem 3.1. In addition, we can conclude that 
the inequalities (7) hold if either is large enough or a is small enough; however, the value of d is not 
presented in (7). In the case that r 2 > d , then small values of K 2 can lead to a unique locally stable interior 
equilibrium. The existences of two or three interior equilibria suggest that under certain conditions, Model 
(6) can have multiple attractors and generate complicated dynamics. 


Theorem 3.3 (Global stability and extinction of one species). Model (6) has global stability at E x 0 = 
(ATi,0) whenever the inequality d > r 2 + holds. Model (6) has global stability at Eqi = ^0, (1 — fq)K 2 ^ 

if '-f < K 2 (1 — ff- J and one of the following conditions is satisfied 


9 

























1. cPnKiKi (h + ^ > rir2(a ^f 1+1)2 ; or 


2 . 

3. 


K\ < 


ah ’ 


or 


Kl > ih and ^hK 1 ^ < 1 ~ if' 


Notes: Theorem 3.3 implies that the necessary condition for the global stability of Eiq is r 2 < d. This 
indicates that the parasite X 2 cannot be facultative if host species X\ persists. If the parasite is facultative 
(i.e., r 2 > d), then the parasite X 2 can drive the host x\ extinct locally when I \2 is large enough and K\ is 
small enough such that the inequalities ifori < tx < K 2 ^1 — hold; local extinction of the host can also 

occur when K 2 is large enough and K±ah > 1 such that the inequality I \2 (l — > ^ 4 ^ 2 holds. 


Ecological dynamics of facultative versus obligate parasites: The condition of the parasite X 2 being 
either facultative or obligate on host Xi has strong effects on the ecological dynamics of the host-parasite 
interaction model ( 6 ). These impacts have been listed in Table 3. We can summarize the main differences 
as follows: 

1. Within the constraints of this model, a facultative parasite is always persistent, because it is not fully 
dependent on the dynamics with its host. If the parasite is obligate, the host always persists. A 
facultative parasite can cause the extinction of the host, while an obligate parasite can go extinction 
itself under certain conditions. 

2. When the parasite is facultative, either large values of or small values of r 2 — d can result in its 
permanence. If the parasite is obligate, both the host and the parasite can persist under intermediate 
values of the parasite death rate d. 

3. When the parasite is facultative, the system is prone to coexistence of both species and can produce 
one, two, or three interior equilibria. When the parasite is obligate, the system can have one or 
three interior equilibria. If the system has two interior equilibria, it exhibits bistability between the 
parasite-only boundary attractor and the interior coexistence attractor. If the system has three interior 
equilibria, it may have two interior coexistence attractors. 


Notes: In the case that the parasite is facultative, i.e., d < r 2 , its population model can be rewritten as 


dx 2 

dt 


= X 2 


eax 1 
1 + hax 1 


+ O 2 - d)(l - 


X2 


w - £) 


which can be rewritten as 


by letting 


dx 2 

dt 


= x 2 


eax 1 X 2 , 

TTteT + r2(1 -i& ) 


V 2 —> r 2 — d and K 2 —> K 2 (l -). 

T'l 

This is equivalent to the case when d = 0 in the host-parasite model. Thus, we can obtain the simplified 
host-parasite model ( 8 ) with Holling-Type II functional responses when the parasite is facultative 


dx 1 
dt 


Xl 



xi \ 

kJ 


ax2 

l-\-hax\ 


dx 2 
dt 


X2 


eax 1 
l-\-haxi 


+ r 2 { 1 



( 8 ) 
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Cases 

facultative parasite ( r 2 > d) 

obligate parasite (r 2 < d) 

Stability of E 00 

Source 

Saddle 

Stability of E w 

Saddle 

Stable if r 2 + x < d 

Stability of f? 0 i 

A saddle if 

^ > K 2 (l — -) > 0: A sink if 

a 2 ^ r 2 J ’ 

^<l-£ 

Does not exist 

Persistence of host x\ 

~cl > K 2 (l A ) > 0 

Always 

Persistence of parasite X 2 

Always 

? ’2 + 1 ZhKi >d>r 2 

Permanence of the system 

n > k 2 (1 - A ) > 0 

a z V r 2 1 

^ + TT^k> d >^ 

Extinction of host x\ 

n > I< 2 (1 - A ) and one of 

the following conditions holds: 

< -k or ( 2 ) A ^i > ik and 

A- 2 (l-^)> ri( 4 1 a + ^ l)2 or 

(°){h 1 d e 'l " r 1 (ahK 1 + l) 2 

r 2 r 2 ) > 3a 2 K±K 2 

Never 

Extinction of parasite X 2 

Never 

d > r 2 + 1+ahK , ■ 

Number of interior 
equilibria 

None, one, two, or three; This 
suggests the possibility of two 
types of multiple attractors: (a) 
two interior attractors; and (b) 
a boundary attractor and an 
interior attractor. 

None, one, or three; This suggests the 
possibility of two interior attractors. 


Table 3: The comparison of sufficient conditions that lead to different dynamical outcomes of a parasite being facultative versus 
obligate for Model (6). 


All the analytic results presented in this section can be applied to Model (8) by letting d = 0. 

Numerical investigations and bifurcation diagrams: According to Theorem 3.2, the ecological model 
(6) can have none, one, two or three interior equilibria depending on parameter values. Unfortunately, the 
explicit forms of these interior equilibria are far too complicated. Our analytical results suggest that the 
values of ?'2 and d have profound effects on the persistence of the system (see Theorem 3.1-3.3) and the 
existence of the number of interior equilibria as well as the stability of these equilibria (see Theorem 3.2). 
For these reasons, we use bifurcation diagrams to explore the effects of r 2 and d on the system’s dynamical 
patterns. 

For convenience, we fix K\ = n = K 2 = 1; h = 4; e = 0.9; a = 5 as a typical example. Figure 1 illustrates 
the number of interior equilibria under this set of parameter values by varying both r2 and d from 0.01 to 
1.5 where the white region means no interior equilibrium; the black region means one interior equilibrium; 
the blue region means two interior equilibria; and the red region means three interior equilibria. The solid 
black line is d = V 2 + below which there is no interior equilibrium. This is supported by Theorem 3.3: 

Model (6) has global stability at Fdo = (Au,0) whenever the inequality d > r2 + 1 + a J^x 1 holds. In Figure 1, 
we can observe that Model (6) has none, one , two, or three interior equilibria when r 2 > d while the system 
has none, one, or three interior equilibria when r 2 < d. This confirms our results in Theorem 3.2 and Table 2. 

Let T 2 = 0.25 (see the dashed black line in Figure 1); we perform an additional bifurcation diagram of the 
stability of interior equilibria with respect to d , as shown in Figure 2, under the same set of parameter values: 
blue means locally asymptotically stable; green means saddle; and red means source. Figure 2 shows that 
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Bifurcation Diagram of d and r 2 



Figure 1: Multiple attractors: Bifurcation diagram of d v.s. r 2 for the number of interior equilibrium when K\ = 7*1 = /C 2 = 
1; h = 4; e = 0.9; a = 5: the while region means no interior equilibrium; the black region means one interior equilibrium; 
the blue region means two interior equilibria; and the red region means three interior equilibria. The solid line is 
d = r 2 + i+ a Jhx 1 = ^2 + 3/14. When the parasite is facultative, i.e., d < T 2 , the ecological system (6) can process none, 
one, two, or three interior equilibria while parasite is obligate, i.e., d > 7 * 2 , the system can process none, one, or three interior 
equilibria. 


when d is small (e.g .,d £ (0,0.1)), there is no interior equilibrium; as d increases (e.g .,d £ (0.1,0.2)), Model 
(6) goes through a saddle node bifurcation (around d = 0.1) which gives one saddle interior equilibrium 
and one stable interior equilibrium along with the local stable boundary equilibrium ^0, K -2 ^1 — ^))’ ^ 

continues to increase (e.g .,d £ (0.2,0.215)), there are three interior equilibria where one is a sink, the second 
is a saddle and the last is also a sink; further increasing d destabilizes the third interior equilibrium (e.g., 
there is a hopf-bifurcation occurring at the third interior equilibrium around d = 0.215); larger d makes the 
system go through a cusp bifurcation (occurring around d = 0.27) which leads to only one stable interior 
equilibrium; at the extreme value of d (i.e., d > r 2 + jT ), the system has global stability at the host-only 
boundary equilibrium (K \, 0). 

One of the interesting questions addressed by the model is what happens to the local stable boundary 
equilibrium ^0, K 2 ^1 — when the system has two interior equilibria after we allow evolution to occur 
(i.e., <Tj > 0,i=l, 2). Can evolution destabilize this boundary equilibrium and save the host from extinction? 
We will explore this analytically and numerically in the next section on the co-evolutionary dynamics of (5). 


4. Co-evolutionary dynamics 


Define 


_ da(ui,u 2 ) 

&U\ — 7T ■ &U2 

OUl 


da(ui,u 2 ) 

d'U2 


Q"U\U\ 


d 2 a(ui,u 2 ) 
0 ^ ’ 
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Bifurcation Diagram on d 



Figure 2: Bifurcation diagrams when K\ = r\ = K 2 = l;r 2 = 0.2; h = 4; e = 0.9; a = 5: blue means locally asymptotically 
stable; green means saddle; and red means source. When the parasite is facultative, i.e., d < 7 * 2 , the ecological system (6) can 
process none, one, two, or three interior equilibria, while the parasite is obligate, i.e., d > r 2 , the system can process none, one, 
or three interior equilibria. In addition, large values of d can destabilize the system and lead to the extinction of the parasite. 


and 


^ U 0 U 0 — 


d 2 a(u\ 1 U 2 ) 


du\ 


Assume that u\,u 2 are trait values such that 


d 2 a(ui,u 2 ) 

du\du2 


d 2 a(u\, U 2 ) 

du\du2 


= K' 2 {u* 2 ) = a Ul {ul,u* 2 ) = a U2 (ut,u* 2 ) = d'(u* 2 ) = 0 


and let 


Ki = Ki(u*), K 2 = K 2 (u 2 ), a = a{u\,u 2 ), and d = d(u 2 ). 


Then according to Proposition 3.1 and Theorem 3.2, the co-evolutionary model (5) can have the boundary 
equilibria (0,0, u*, u 2 ), E Xl0uiU2 = (K 1 , 0, n*, u 2 ), E 0x2UlU2 = (o,K 2 (l - and potentially mul¬ 

tiple interior equilibria (xl,x 2 ,ul,u 2 ) depending on the values of K i: ri,i = 1,2 and d,a. In general, the 
following theorem provides the existence and local stability results on the facultative-parasite-only equilib¬ 
rium Eq X 2UiU2 and the host-only equilibrium E Xl q U iU2 for the co-evolutionary model (5): 


Theorem 4.1. The host-only equilibrium E Xl q U iU2 = (Ki,0,ul,u 2 ) exists if 

v , ( n j iir eKr{ul)a U2 {ul,ul) 

R 1 (u 1 ) = 0 and d (u 2 ) = 

(1 + ha{u l ,u 2 )Ki{u 1 )Y 

and it is locally asymptotically stable if the following inequalities hold 

ea(u\,u 2 )K\(u\) 


0 < 


1 + ha(u \, u 2 )Ki(u\) 


<d(u*)-r 2 , K"(u\) < 0, 
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and 


eKxjup [a U2U2 {ul,ul)(l + haKijul)) - 2Ivi(u*)<4 2 {u{, u 2 )\ 
(1 + ha(u \, U 2 )/ii(it*)) 3 

The parasite-only equilibrium Eq X2U1U2 = ^ 0,/\ 2 ^1 — exists if 

r 2 K' 2 (u* 2 ) 
I< 2 (u* 2 ) 


< d"{u 2 ). 


a Ul (ul,u 2 ) = 0, d{u 2 ) < r 2 and d'(u 2 ) = 


and it is locally asymptotically stable if the following inequalities hold 


r i 


a(ul,u* 2 ) 


<K, 


! (« 2 *)(l-^), a UlUl {u\,u* 2 )> 0, 


and 


r 2 (l - ^) [K 2 {u* 2 )K'f{u* 2 ) - 2{K' 2 {u* 2 )Y 


< d'\u* 2 ). 


Notes: According to Proposition 3.1, the extinction equilibrium Eqq is alway unstable for the ecological 
model (6). The host-only equilibrium E 10 is locally asymptotically stable when r 2 + < This 

condition implies that social parasite should be obligate, i.e., r 2 < d. The facultative-social-parasite-only 
equilibrium Eoi is locally asymptotically stable when the inequality < 1 — ^ holds. The results of 
Theorem 4.1 imply that the local stability of the facultative-social-parasite-only equilibrium Eq X 2UiU2 and 
the host-only equilibrium E Xl o UlU2 are determined by the concavity of the trait function Ki(vi ), i = 1,2 and 
a(vi,v 2 ) evaluated at the equilibrium under conditions that the equilibrium is Ecologically Stable (i.e., ES). 
If the death rate due to parasitism d is independent of the trait value u 2 , then we have the following corollary: 


Corollary 4.1. The host-only equilibrium E Xl o UlU2 exists if A'((u*) = a U2 (ul,u 2 ) = 0, and it is locally 
asymptotically stable if the following inequalities hold 


0 < 


ea(u\, u 2 )K\(u\) 

1 + ha(u \, u 2 )Ki{u\\ 


<d — r 2 , K"(ui)<0, and a U2U2 (ul,u 2 ) < 0. 


The parasite-only equilibrium Eq X 2UiU2 exists if K 2 (u 2 ) = a Ul (ul,u 2 ) 
stable if the following inequalities hold 


< I\ 2 (u* 2 ) (l- — 


Q'U 1 ui('U'h U 2 ) > 0 , 


= 0, and it is locally asymptotically 

and K 2 (u 2 ) < 0. 


f.l. Boundary equilibria and ESS 

The trait functions can have many forms, such as Gaussian distributions, polynomial, exponential func¬ 
tions (Abrams 1990; Bergelson et al 2001; Mostowy and Engelstadter 2011; Nuismer et al 2012; Landi et 
al 2013). In this subsection, we apply the results of Theorem 4.1 and its Corollary 4.1 to some specific 
trait functions to explore how these functions affect whether or not the boundary equilibrium E Xl q U iU2 or 
Eox 2 u 1 u 2 can have ESS. More specifically, we assume that 

— Vi — U 2 + c ) 2 ( V 1 ~ v 2 — c ) 2 (*> 2 - c ) 2 ( y 2+°) 2 

Ki(vi) = K 0i e 2 ° K i , a(vi,v 2 ) = a 0 e , and d(v 2 ) = d 0 e ^ ,*=1,2 

where c > 0 and Vi € [0, c]. These chosen trait functions are modified from Gaussian distributions. 
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4-1.1. The fixed parasite death rate 


In this subsection, we apply the results of Corollary 4.1 to particular trait functions when the death rate 
d of the parasite due to attacking all potential hosts is independent of the trait U 2 - Let trait values set be 

Uj = [0, c], * = 1, 2. We assume that Ki(vi) = K 0i e 2 " K i ,i = 1,2 and a(v i, u 2 ) = age 
which gives follows: 


( U 1 -w 2 +c) 2 (^i _^ 2 _ c ) 2 

XT - 


and 


_ 

A”'0) = Ki(c) = Kqi , a(0,0) = a(c, c) = age ^, a(0,c) = a(c, 0) = ag,i = 1,2 
/v'(0) = K[ (c) = a Vi (0, 0) = a Vi (c, c) = a Vi (0, c) = a Vi (c, 0), * = 1,2 

_ c 2 

A'"(0) = K"(c) = < 0 , a ViVi (0,0) = a ViVi (c,c ) = > 0, 

a Ki °a 

O'ViVi ( 0 ? c) — dviVi (C 0 ) 9 ^ 0 . 

at 


This implies that the trait dynamics of the co-evolutionary model (5) are positively invariant in the trait 
space Ui x U 2 = [0,c] x [0, c]. 


The parasite-only equilibrium: According to Corollary 4.1, the parasite-only equilibrium Eg X2UlU2 = 
(0, Kg 2 (1 — ) , ul, u 2 * ) is locally asymptotically stable (i.e., CS) if the following conditions hold 


K' 2 {ul) = a ui (ul,u 2 ) = 0, 


r 1 


1 - 


< aoKo 2 e 


U1 (ut,u2) = 


aoc 2 e 


> 0, and K 2 (u 2 ) < 0 


( 9 ) 


where (ul,u 2 ) = (0,0) or (ul,u 2 ) = (c, c). Then the fitness functions of host and parasite at Eq X2U1U2 are 
represented as follows: 


Gi(ui,u*,x*) = r 1 -a(v l ,u* 2 )Kg 2 [l-=r 1 -agKg 2 {l--^je 

0 -*) 


G 2 {v 2 ,u*,x*) = —d + r 2 (l — 


Kq2(i-^ 

K 2 (v 2 ) 


(vi - u % + c) 2 (« 1 - c) 2 


„ 2 („ 2 -c) 2 


= —d + r 2 1 — (1 — 


which gives 


max {Gi(ui, u*, £*)} = Gi(it*, u*, x*) 

uiGUi 


ri - a 0 K 02 



< 0 


and 


max {G 2 (i> 2 , «*,£*)} = G 2 (u 2 ,u* ,x*) 

v 2 £U 2 


-d + r 2 




= 0. 


Therefore, according to the definition of ESS (i.e., (2)), we can conclude that the two parasite-only equilib¬ 
rium ^0, A 02 (l — ,0,0^ and ^0 ,/\q 2 ^1 — , c, of the co-evolutionary model (5) are ESS. 


Biological scenarios: The case studied above suggest that the parasite-only-equilibrium Eg X2UlU2 can have 
two ESS strategies (u*,u 2 ) = (0,0) and (■ u\,u 2 ) = (c, c) when the parasite is facultative (i.e., d < r 2 ). This 
may describe the case of facultative slavemakers, such as Formica subnuda whose colonies can survive as 
slaveless (Savolainen and Deslippe 1996). 
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The host-only equilibrium: According to Corollary 4.1, the host-only equilibrium E Xl q U iU2 = (A'oi, 0, u\, u 2 ) 
is locally asymptotically stable if the following conditions hold 


K[(u\) = a U2 (ul,U 2 ) = 0, 


eaoKoi 
1 + haoKoi 


< d — T 2 , a U2 u 2 (ul,u 2 ) = — 


2aoc 2 


< 0, and K '2 (^ 2 ) < 0 


( 10 ) 


where {u\,u 2 ) = (0,c) or {u\,u 2 ) = (c, 0). Then the fitness functions of host and parasite at E Xl q U iU2 are 
represented as follows: 


Gi{vx,u*,x*) = n 1 - 


G 2 (v 2 ,u*,x*) = 


K 01 


K 0 ie K i 


= n 1 - e 


v ~\ 01“ c ) 

72 — 


ea{u\,V2)Ko 1 


1 +ha (ul,v 2 )K 01 


— d + r 2 = — d + r 2 + 


ei^oiape 


(^T -u 2 +c)^(n| -v 2 -c) z 


l+ZiKoiape 


(^^-u 2 +c) 2 (^*-u 2 -c )2 

2 rr 4 


which gives 


max {Gi(ui, u*, x*)} = Gi(rq,u*,x*) = 0 

uieUi 


and 


max{G 2 (u 2 ,M*,x*)} = G 2 (u 2 ,u*, 

v 2 eu 2 


**) = 


d + r 2 + 


ea 0 A 01 

1 + hag Agi 


< 0. 


Therefore, according to (2), the two host-only equilibria (A' 0 i, 0, 0, c) and (A 0 i, 0, c, 0) of the co-evolutionary 
model (5) are ESS. 


Biological scenarios: The cases studied above suggest that for the trait functions given, the host-only equi¬ 
librium E Xl o UlU2 of the co-evolutionary model (5) can have two ESS (w*,rt 2 ) = (0,c) and ( u*,u 2 ) = (c, 0) 
when the parasite is obligate (i.e., d > r 2 ). This can be classified as one co-evolutionary outcome when the 
host successfully resists invasion by the parasite, resistance occurs via effective front-line defenses (Kilner 
and Langmore 2011). For example, Ortolani and Cervo (2010) show that Polistes dominulus foundresses are 
now so large and aggressive that they consistently defend their nests from attack by the brood parasite P. 
sulcifer, and that parasitism is rarely seen. 


4.1.2. The parasite death rate depending on its trait 

Our application of Corollary 4.1 in the previous subsection indicates that both boundary equilibria 
E X1 Ouiu 2 an d Eq X2U1U2 cannot simultaneously have ESS strategies, due to the fact that d and r 2 are fixed. 
The interesting questions become, if the death rate d of the parasite also depends on its trait u 2 , can Model 
(5) have both boundary equilibria E Xl o UlU2 and Eq X 2UiU2 being locally asymptotically stable (i.e., 
CS) under different trait values of u 2 ? If it can, is it possible for these two boundary equilibria 
to have ESS? 

To investigate these questions, we should apply the results of Theorem 4.1 by letting d depend on the trait 

_ (-»2— c ) 2 (^2+ c ) 2 

value u 2 . To continue our study, we choose d{v 2 ) = d^e 2,7 ^ which has the following properties: 

c 2 

d( 0) = doe^d d(c) = do, d'( 0) = d'(c) = 0, d"(0) = ~ J ^° r \ —— > 0, and d"(c) = — < 0 . 

a d a d 

If we assume that Ki(yi),i = 1,2 and a(vi,v 2 ) have the same trait functions as before, then we can conclude 
that the trait dynamics of the co-evolutionary model (5) are positively invariant in the trait space HJi x U 2 = 
[0,c] x [0,c]; and the equilibrium trait values (u\,u 2 ) are u* = 0 or c for both i = 1,2. According to 
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Theorem 4.1, we can conclude that both boundary equilibria E Xl o Ul u 2 = (Aoi,0, 0,c) and Eq X 2 Ui u 2 = 
( 0 , Kq 2 (1 — 77 ) 10 , 0 j have convergence stability if the following conditions hold: 


d 0 <ri(l + a 0 hK 0 i) ea 0 K 01 


< 


1 + ag/iAoi 


< do — T 2 and r± < aoA'o 2 e 



For example, if we let 


ri = 0.1, r 2 = 1.5 < do = 2.1, a 0 = 2, e = .9, h = c = a k 1 = cfk 2 = A '01 = 1, A' 02 = 100, ad = 1.05, a a = -56, 

then both E Xl o UlU2 = (A' O i,0, 0,c) and E 0 x 2 u 1 u2 = (o, A" 02 (l — yj ) , 0,c) are locally asymptotically stable 
(CS). However, the equilibrium trait (zt*,u 2 ) = (0,c) is not ESS for E Xl o Ul u 2 = (AToi, 0,0, c); and (u*,u 2 ) = 

(0,0) is not ESS for E 0x2UlU2 = ^0, K 0 2 ^1 — d ° e r2 d ^ , 0, 0^ . See the fitness functions Gi(vi,u*,x*),i = 

1,2 for these two boundary equilibria shown in Figure 3-4: Figure 3(a) and 4(a) show that the strate¬ 
gies (ui,u 2 ) = (0,c) and (uf,^) = (0,0) are ESS for the host at the boundary equilibrium (AT O i,0), 


0, Ki 


02 


1 - 


dpe 


respectively. However, Figure 3(b) and 4(b) show that these strategies are not 


ESS for the parasite at these boundary equilibria, since there are trait values v 2 such that 62 (^ 2 , u*, x*) > 0. 




(a) The fitness function of host (i.e., Gi(vi,u* ,x*)) with u* = (b) The fitness function of parasite (i.e., G 2 {v 2 ,u* ,x*)) with 
(0, c) and x* = (ifoi,0) u* = (0, c) and x* = (/Coi,0) 

Figure 3: The fitness functions of host (i.e., Gi(vi,u*, x*)) and of parasite (i.e., G 2 {v 2 ,u*,x*)) for the co-evolutionary model (5) 
when n = 0.1, 7*2 = 1.5 < do = 2.1, ao = 2, e = 0.9, h = c = aK 1 = <xk 2 — -^01 = 1, K 02 = 100, ad = 1.05, a a = 0.56. The 
strategy u* = (0,0) is not an ESS, because there are trait values V 2 such that G 2 (v 2 ,u*, x*) > 0. However, the co-evolutionary 
system has local stability at ( u*,x *). 
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Figure 4: The fitness functions of host (i.e., Gi(vi,u*, x*)) and of parasite (i.e., G 2 (u 2 , u *, x*)) for the co-evolutionary model 
(5) when ri = 0.1, r -2 = 1.5 < do = 2.1, ao = 2, e = 0.9, h = c = ctk 1 = &k 2 = tToi = 1, K 02 = 100, <?d = 1.05, cr a = 0.56. 
The strategy u* = (0, c) is not ESS for since since there are trait values V 2 such that (J 2 {y 2 , u *, x *) > 0. However, the 
co-evolutionary system has local stability at ( u*,x *). 


Biological scenarios and implications: The case studied above suggests that a host and parasite can 
employ multiple strategies that generate convergence stability but not ESS. This may be due to the in¬ 
teractions of different strategies and counter-strategies acquired at different stages of co-evolution (Kilner 
and Langmore 2011). For example, parasites may initially be facultative, then acquire host-specific signa¬ 
tures, but through counter-selection by hosts may subsequently revert to genetic signatures expressed before 
parasitism. As one possible example, among the insects, selection to become chemically insignificant may 
have facilitated the ability to acquire host hydrocarbon signatures previously after parasitism (Kilner and 
Langmore 2011). Such a case could be applied to slave-making ants, which can be either obligate social par¬ 
asites, depending on enslaved hosts ants throughout their whole lives (Topoff and Zimmerli 1991; Ruano et 
al. 2013) or alternatively facultative slave-makers. Facultative slave-making ants, like those in the Formica 
sanguinea complex, may engage in slave making, but individual colonies are able to revert to producing 
their own workers if parasitized workers are removed (Topoff and Zimmerli 1991). They could represent an 
intermediate parasitic group, between freeliving species on the one hand, and obligate slave-making species 
on the other. 


4-2. Trait functions follow Gaussian distributions 

To continue our analytical study, we focus on the co-evolutionary dynamics of Model (5) for chosen trait 
functions of Ki{vi) and a(v 1 , 1 x 2 )- For the rest of the section, we let the death rate d of the parasite X 2 be 
independent of the trait 1 x 2 ; and we assume the following functional forms for the carrying capacity Kiivi) 


18 
















and the parasitism efficiency a(v \, v 2 ): 


Ki{vi) 

= K 0i e 

a(v i,v 2 ) 

(»i -«a) 

= a 0 e 


( 11 ) 


a(v!,v 2 ) denotes the parasitism efficiency of a parasite with phenotypic trait v 2 on host individuals with 
phenotypic trait v\ with the assumption that the stronger host-parasite interactions are, the more similar 
host and parasite traits are. The symmetric form of a in (11) has been previously used in the study of 
character displacement by Taper and Case (1992), which assumes that the parasitism efficiency is normally 
distributed around a maximum value of ao with a variance as a function of trait difference in the parasite 
and host vi — v 2 . The larger value of cr^, the greater sensitivity of the parasitism efficiency a changes with 
respect to the changes in trait difference vi — v 2 . There are other alternative functions of a (see Dieckmann 
and Marrow 1995; Doebeli and Dieckmann 2000; Zu et al 2007), for example, the asymmetrical predation 
efficiency a which has been previously used in the study of character displacement can be described as 


a{v i,v 2 ) 


a 0 e~ 


--- ^ 2rr 2 


( 12 ) 


where (3^0 (Doebeli and U. Dieckmann 2000). We also assume that the resource availability for the host 
and parasite varies with their phenotypic trait Vi, which follows a Gaussian distribution iV^O,^.). This 
distribution assumes that the maximum inherent equilibrium level Kqi for each individual single species Xi , 
using strategy Vi, in the absence of other species, is attained at trait Vi = 0, and that Ki{vi) is normally 
distributed around Vi = 0 with variance af c .. The larger variance crj^., the greater sensitivity of the carrying 
capacity Ki(vi) changes with respect to the changes in the trait Vi. Fix the trait value of the parasite u 2 , 
the formulations of K i(vi) and a(v\,u 2 ) indicate that increases in the host’s trait value |ui| result in the 
decreased values of K\(v{) and a(vi,u 2 ) when v\ < 0 or V\ > u 2 but can result in the increased values of a 
when i>i € ( 0,u 2 ). This implies that there is a trade-off between the carrying capacity of the host and the 
parasitism efficiency for a certain range of v\ when v\ < 0 or v\ > u 2 . On the other hand, for a fixed host’s 
trait value u\, there is no such trade-off for the parasite when v 2 < 0 or v 2 > u\ but a trade-off exists when 
v 2 e (0,ui). 


As a consequence, the co-evolutionary dynamics of the monomorphic resident host and parasite popula¬ 
tions with residence traits of u\ and u 2 in the host, parasite, respectively, are given by the following set of 
nonlinear equations: 


dx 1 
dt 


x\Gi{vi,u, a;)|„ 1=ui = xiH 1 (u 1 ,u 2 ,x 1 ,x 2 ) 


Xi 



Xl \ 
Kl(ui)/ 


a(u 1 ,u 2 )x 2 
l-\-ha(u\ ,u 2 )xi 


dx 2 
dt 


dui 

dt 


x 2 G 2 (v 2 ,u , x)\ V 2 —u 2 — x 2 H 2 (u 1 ,u 2 ,x 1 ,x 2 ) 


x 2 


ea(u\ ,u 2 )x\ 
l-\-ha(ui ,u 2 )xi 


d + r 2 (l 


%2 ) 
K 2 (u 2 )' 


„2 dGi(vi,u,x) | _ „2 dHi (ui,u 2 ,xi,x 2 ) _ 2 

1 dvi \vi=ui 1 dui 1 


; • i x i u i 

i(“i) 


+ 


(ui—u 2 )x 2 a(ui,u 2 ) 
a‘a(l-\-ha(ui ,U 2 )xi ) 2 


(13) 


du 2 

dt 


2 dG 2 (v 2 ,u,x) I _ 2 dH 2 {u 1 ,u 2 ,xi,x 2 ) 2 

2 dv 2 \v2=u 2 2 Qu 2 2 


r 2 .x 2 .u 2 1 e(m— u 2 )xia(ui,u 2 ) 

!?'k 2 K 2 (u 2 ) a-Z(l+ha(u 1 ,u 2 )x 1 ) 2 ’ 


where we take a with the symmetric form given in (11). Since the death rate of the parasite x 2 does not 
change over time, thus the difference of r 2 and d determines whether the parasite is facultative or obligate: 
the parasite is facultative when r 2 > d > 0, and is obligate when d > r 2 > 0. Our model (13) allows us to 
investigate the ecological and evolutionary dynamics of the host and parasite interactions, where the parasite 
can be either obligate or facultative. In particular, Model (13) includes the special case which is studied by 
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Zu et al. (2007). 


Assume that (xj, x 2 , u^, u 2 ) is an equilibrium of the co-evolutionary host-parasite model (13). We say 
(x*, x 2 , u*, u 2 ) is a boundary equilibrium if x\x 2 = 0 while it is an interior equilibrium if x*X 2 > 0. The 
equilibrium (x*, x 2 , u*, u 2 ) satisfies the following equations 


dx i x i \ a(iti,u 2 )x 2 

—— = 0 => xi = 0 or n 1 — ———- = --— t- -— 

at \ Ai(ui)/ 1 + na{ui, u 2 )xi 

„ „ ea(ui,u 2 )xi / x 2 \ 

_ = „^ X2 = 0or 1 + Muiitt2)xi+ r^l--^- y j= (i 

dui _ riXiUi _ (iti - u 2 )x 2 a(u 1 ,u 2 ) 

dt aj < - i Ki(ui) a 2 (l + ha{ui,u 2 )xi) 2 ’ 

rfu 2 _ r 2 x 2 u 2 _ e{u\ — u 2 )x\a(ui,u 2 ) 

dt <7 2 K2 K 2 {u 2 ) 0-2(1 + ha(ui,u 2 )xi ) 2 ' 


(14) 

(15) 

(16) 
(17) 


We have the following proposition regarding the boundary equilibria of (13) and their stability. 


Theorem 4.2. [Equilibi'ia of the evolutionary model (13)/ The evolutionary model (13) always has 
the following two boundary equilibria Eoooo = (0,0,0,0), E Xl ooo = (Aoi, 0,0,0) where E Xl ooo is locally 
asymptotically stable if r 2 + i+°hK 01 < d and it is a saddle if r 2 + > d. In addition, the following 

statements are true: 

1. If r 2 > d, then Model (13) has the third boundary equilibrium Eo X2 qq = (0,Kq 2 ^1 — ,0,0). Both 

Eoooo and Eo X2 oo are always saddles. 

2. Model (13) can have one, two, or three interior equilibria (xj, x 2 , 0,0) where (x^x^) is an interior 
equilibrium of the ecological model (6) by letting K\ = Kqi, K 2 = Kq 2 and a = ag. 


3. The interior equilibrium (x*,X2,0, 0) is locally asymptotically stable if the following conditions hold: 


m { K 02 ’ v'ic Koi } > a °- a- 0 }a ? 02 + 


max 


2 r 2 

K 02 

/ °~ K 2 


(1+ao/iiCoi) 3 


> alhK Q2 (r 2 + f), and 


n eKn 2 
T2 Koi ' 


J K 2 


aoKo 


—} > % 
l J 


(18) 


Notes: According to Corollary 4.1, one necessary condition for Eo X2Ul u 2 being locally asymptotically sta¬ 
ble is that a UlUl (it*, u 2 ) > 0. The results of Theorem 4.2 imply that the boundary equilibrium Eo X2 oo = 
(0, AT 02 ^1 — ,0,0) cannot be locally stable, since a„ iui (0,0) < 0. Considering the result in Theorem 

4.1 indicates that the trait function of parasitism efficiency a(ui,u 2 ) plays an important role in determin¬ 
ing whether the boundary equilibrium Eq X2 oo can be locally stable or not. In addition, Condition (18) in 
Theorem 4.2 indicates that a large ratio of the host (or parasite) intrinsic growth rate r± to the maximum 
carrying capacity of host (or parasite ) A'oi, a large ratio of variance of the carrying capacity of parasite 
to host, i.e. and large values of the variance of the trait difference of parasitism efficiency a, can lead 
to local stability of the coexistence of host and parasite (xj, x 2 , 0,0). According to Theorem 3.1, for fixed 
values of n,Koi, a 0 , e, the proper value of d can guarantee the permanence of the ecological system, and thus 
guarantee the existence of the interior equilibrium (x*,x 2 ). In the rest of this subsection, we will focus on 
the ESS of the host-only boundary equilibrium E Xlt 0 , 0,0 and the interior equilibrium (xJ,X2,0,0). 


Theorem 4.3. [ESS of co-evolutionary model (13)/ The strategy (u\,u 2 ) = (0,0) of the boundary equi¬ 
librium E Xl 000 is an ESS if the following inequality holds 


r 2 


eaoAoi 
1 + aohKoi 


< d. 
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The strategy (itj/u^) = (0,0) of the interior equilibrium (a^, £ 2 ? 0,0) is an ESS if Condition (18) holds and 
the following inequalities hold 


r i 

< J K 1 K 01 


> 


max 


I a K 2 j r 1 eK 02 ' 

\ a K 1 V r 2 K 01 ’ 


7 k 2 ° 


(19) 


Notes: By simple rearrangements, we have the following two equivalent equations: 


r i 


^h K or 


> 


<?K-, 


£l eK m 

r 2 K 0 i 


a 0 




o'/fi c k 2 


> 


eK(ji K[y2 

nr 2 


-do 


and 


r i 




d^A'oi 


> 


apKoi 


-ao 


o- Kl cr K2 < 1. 


Thus, according to the proof of Theorem 4.2 and 4.3, we have the following corollary: 
Corollary 4.2. Assume the following inequalities hold 


. f 2 r 2 ri j 

m <-, —^-> 

l Kq 2 & is- Km j 


<7 2 Kl K 01- 


> a o, 


r\ r 2 


Kq\Kq 2 (1 + a 0 hK 01 ) 3 


> alhK 02 (r 2 + -)• 


The strategy = (0,0) of the interior equilibrium (x*, x%, 0,0) is an ESS if one of the following 

conditions holds 


i. 1 * 2 . /negsz > 4 and > J 

<?K 1 V r 2*M)l X 1 VK 1 <TK 2 V 

and <Jk 1 o’k 2 < 1 - 


eK oi K 02 

r 1 r 2 


a o- 


2 . 


r l°K 2 ° 

CLqKqi 


> 


Notes: By comparing the results of Theorem 4.2 with Theorem 4.3 and its corollary 4.2, we can conclude 
that the variances of the proposed trait functions, i.e., 07^,0 t k 2 and o^, play essential roles in guarantee¬ 
ing the local stable interior equilibrium (a;*, x 2 , 0, 0) having an ESS. More specifically, large variances of 
the parasite carrying capacity <tk 2 and of the parasitism efficiency cr Q are required to make sure that the 
strategy {u\,u 2 ) = (0,0) of (x*, x%, 0, 0) is an ESS but it may not be global ESS because it is possible that 
the system has other ESS strategies for coexistence. This can be considered as a tolerance of the parasite 
co-evolutionary outcome, where hosts not only concede to the parasite and accept it in their nests, but also 
make adjustments to their life history (or other traits) to minimize the negative effects of parasitism on their 
fitness (Kilner and Langmore 2011). This co-evolutionary outcome is more likely be the case that complete 
parasitic control of the co-evolutionary trajectory. This could fit with the ecological findings of imperfect 
recognition of parasite eggs, and a low but positive level of acceptance of brood parasitism by some cuckoo 
and cowbird host species (Davies et al 1996; reviewed by Kruuger 2007). 


Theorem 4.4. [The unique ESS of the co-evolutionary model (13)/ The strategy (u\,u 2 ) = (0,0) of the 
boundary equilibrium E Xl 000 is the unique ESS of the co-evolutionary host-parasite model (13) if the following 
inequality holds 


eap-Koi 

1 + ahKoi 


< d. 


Assume that <jk 2 > then the strategy {u\,u 2 ) = (0,0) of the interior equilibrium {x\, 4,0,0) is the 
unique ESS of the co-evolutionary host-parasite model (13) whenever the following conditions hold 


Kq 2 O .0 < 1 KoiOq < ^ 

■> „ n r l r 2 _I_ ea O _ 

O’ Kq\Kq 2 (1+ao/iiCoi) 3 

aK i j rfeKpf £2 °a \ j eKgiKgf 
0-ftTi y r 2 K 01 > a Kl a K2 \J r x r 2 


f 2 r 2 n 1 

\ K<J 2 ’ tr'x K 01 l 


> alhK Q2 {r 2 + f) and 
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Notes: Theorem 4.4 provides sufficient conditions when Model (13) has a unique ESS. One direct impli¬ 
cation is that if (it*,«;)) = (0,0) is an ESS of the host-only boundary equilibrium -E^ooo, then Model (13) 
cannot have other locally asymptotically stable equilibrium. Thus a large death rate due to the parasite 
hunting/attacking the host can lead to the extinction of the parasite, and lead to a successful resistance by 
hosts evolutionary outcome (Kilner and Langmore 2011). 

On the other hand, when the ratio of the variance of the parasite to the host is larger than 1, i.e., > 1, 

small values of a 0 can lead to the interior equilibrium (xj, x 2 , 0 , 0 ) being the only equilibrium of the co¬ 
evolutionary model (13) having ESS. This type co-evolutionary outcome is referred to as acceptance of the 
parasite , which can be considered an adaptive strategy for hosts when the costs of rearing a parasite are, on 
average, lower than any recognition costs (reviewed by Kruuger 2007; Kilner and Langmore 2011). 


Our analysis suggests that it is possible for Model (13) to have multiple interior equilibria with potential 
ESS in the following scenarios: 

1. The condition > gk 1 does not hold in Theorem 4.4. 

2. The ecological model can have multiple stable equilibria which also can have ESS. 


We have been focusing on the equilibrium (xf, x 2 , u 2 ) °f the co-evolutionary model (13) when u\ = 

u 2 = 0- It is possible that our model (13) has an equilibrium with u\u 2 7 ^ 0. See the following proposition: 

Proposition 4.1. If (x*, x 2 , u^, u 2 ) is an equilibrium of the co-evolutionary model (13) with u\u 2 7 ^ 0, then 
we have x\x 2 > 0 , and x|,x 2 , Wi,m 2 can be solved from the following equations: 


where 


x *i = fi{u*,u 2 ) = 

g(ui,u 2 ) 

fl(Ul,U 2 ) 

/ 2 (mi,m 2 ) 


g(u\,u* 2 ) 


= Mu*,u 2 ) and 


r\u l 


_ = K ~ u 2 )g(ul,u 2 )a{ut,u 2 ) 

4/lK) 0-2 (1 + Mu*,W*)/i(-uJ,W*)) 2 


eriff| 2 tiiK2(ti2) 
r 2°K 1 “2.K'l(ui) 

Ki _ ry+aKig , \/_(r 1 +aK 1 g) 2 +r 1 ahK 1 (r 1 ahK 1 +2r 1 -2K 1 g) 

2 2 r±ah _ 2 /'i ah 

K 2 j_ aeK2 
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2r2dhg 


1 yj (r2g—aeK2) 2 +r2ahK2{r2ahK2+2r2g+‘2eaK2) 

2 ah 2 r\ah 


Notes: Since our trait functions are even function, so are fi(u\,u 2 ) and g(u\,u 2 ), i.e., 

fi(u 1 ,u 2 ) = f 1 {-u 1 ,-u 2 ),f 2 (u 1 ,u 2 ) = / 2 (-u 1 ,-u 2 ) and glu^u^) = 5 (-m 1 ,-u 2 )- 


Therefore, if {u\,u 2 ,x\,x 2 ) is an interior equilibrium of Model (13), then so is (—— u 2l x*, x 2 ). Propo¬ 
sition (4.1) suggests that the co-evolutionary model (13) could have multiple interior equilibria, and have 
potential multiple interior ESS, with different trait values. Since we are not able to solve the interior equi¬ 
librium (x*, x 2 , ul, u 2 ) explicitly for the case of u\u 2 7 ^ 0 , we use numerical simulations to illustrate typical 
scenarios of the interior equilibrium (x^,x 2 ,wj, u 2 ) when u\u 2 7 ^ 0. We are particularly interested in the 
co-evolutionary outcomes when the parasite is facultative, i.e., r 2 > d. As an example, we fix the values of 
parameters as follows 


Kqi = rq = Kq 2 = 1; 00 = 5; h = 4; e = 0.9; r 2 = 0.25 > d = 0.185. 

We would like to note that the chosen parameter values represent typical dynamics of Model (13) when the 
ecological model ( 6 ) processes two interior equilibria, e.g., the blue region of Figure 1. 
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According to Figure 2, we can see that, under this set of parameter values, the ecological model (6) has one 
saddle interior equilibrium; one stable interior equilibrium (0.5428,1.0841); and one stable boundary equilib¬ 
rium (0, 0.26). Depending on initial conditions, the trajectory of (6) will converge to either (0.5428,1.0841) 
or (0,0.26). For example, if we take (aq(0), £ 2 ( 0 )) = (0.5,2) as the initial conditions, the ecological model 
(6) converges to the parasite-only boundary equilibrium (0,0.26). After we turn on evolution, we have the 
co-evolutionary model (13). According to Corollary 4.1 and Theorem 4.2, the parasite-only boundary equi¬ 
librium (0,0.26,0,0) is unstable. 

To investigate the co-evolutionary outcome, we fix the initial conditions (aq(0), £ 2 ( 0 ), ui(0), U 2 (0)) = 
(0.5,2,1,0.1), let cr_R' 1 = <jk 2 = 1, and vary the value of <r a . More specifically, we simulate the dynamics of 
the co-evolutionary model (13) when cr a = 2, 1, 0.5, 0.15, 0.05, 0.005. The simulations are shown in Figure 
5-9 where the host X\ is red; the parasite X2 is blue; the trait value for host u± is green; and the trait value for 
parasite U 2 is cyan: When er Q = 2, the co-evolutionary model (5) converges to the stable interior equilibrium 
(0.0102, 0.3452, 2.3846,0.1287) locally (see Figure 5); When a a = 1, the co-evolutionary model (5) converges 
to the stable interior equilibrium (0.0619,0.4984,1.8290,0.4813) locally (see Figure 6); When a a = 0.5, the 
co-evolutionary model (5) converges to the stable interior equilibrium (0.5428,1.0841,0,0) (see Figure 7). 
However, when cr a = 0.15, 0.05, the co-evolutionary model (5) has fluctuating dynamics (see Figure 8-9). 
And when the value of cr a is small enough, e.g., o a = 0.005, the co-evolutionary model (5) has fluctuating 
dynamics first but eventually converges to the stable interior equilibrium (0.5428,1.0841,0,0) (see Figure 
10). These simulations suggest follows: 



t 


Figure 5: The stable interior equilibrium is (0.0102,0.3452,2.3846,0.1287). The dynamics of the co-evolutionary model (5) 
when uki = (J k 2 ~ Koi = n = ^02 = 1, ^2 = 0.25 > d = 0.185, ao = 5, e = 0.9, h = 4, a a = 2 with initial values 
(iri(0),#2(0),ui(0), 1^2(0)) = (0.5,2,1,0.1): the host x\ is red; the parasite X 2 is blue; the trait value for the host u\ is green; 
and the trait value for the parasite U 2 is cyan. 


1. The stable interior equilibrium (xi,X2,ul,u%) with u*u* 2 ^ 0: The co-evolutionary model (13) 
can indeed have locally stable interior equilibrium (aff, x|, ul, u%) with 7 ^ 0 (see Figure 5 and 
6 ). This occurs when the ecological dynamics ( 6 ) has two attractors in the absence of evolution: the 
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Figure 6: The stable interior equilibrium is (0.0619,0.4984,1.8290,0.4813). The dynamics of the co-evolutionary model (5) 
when ctkx = °’k 2 ~ i^oi = r\ = Kq2 = 1, r2 = 0.25 > d = 0.185, ao = 5, e = 0.9, h = 4, a a = 1 with initial values 
(rci(0), ^ 2 ( 0 ), ui(0), 7 x 2 ( 0 )) = (0.5,2,1,0.1): the host x\ is red; the parasite X2 is blue; the trait value for host u\ is green; and 
the trait value for the parasite U 2 is cyan. 

parasite-only equilibrium (o, A '02 (l — yj) ) and the interior attractor where both the host and parasite 
coexist. 

2. ESS or not: Calculations show that: (a) when a a = 2, 1, the co-evolutionary model (5) converges 
to the stable interior equilibrium (0.0102, 0.3452, 2.3846,0.1287), (0.0619, 0.4984,1.8290, 0.4813) locally, 
respectively, which are not ESS; and (b) when a a = 0.5, 0.005, the co-evolutionary model (5) converges 
to the stable interior equilibrium (0.5428,1.0841, 0,0) where the strategy (0, 0) is an ESS. 

3. Evolution can save the host from local extinction by creating new stable interior equilibria 
(see Figure 5-6) or shifting the basins of attractions (see Figure 7): In the absence of evolution, 
under the condition of A' 0 i = rq = A '02 = 1; ao = 5; h = 4; e = 0.9; r 2 = 0.25 > d = 0.185 and 
(aq(0),£ 2 ( 0 )) = (0.5,2), the ecological model (6) converges to the parasite-only equilibrium (0,0.26); 
However, in the presence of evolution, the co-evolutionary model (5) can converge to the stable interior 
equilibrium locally where both the host and parasite can coexist. 

4. Effects of the variance er a of the trait difference in the parasitism efficiency a(ui,U2): Small 
values of a a can destabilize the evolutionary dynamics and generate oscillations between population and 
trait dynamics (see Figure 8-9), while extremely small values of cr a can cause catastrophe events where 
the large oscillations collapse with the boundary of the basins attractions of the stable interior equi¬ 
librium (:r*, X 2 , 0,0) with the consequence that the trajectory converges to {x*, x%, 0,0) (see Figure 10). 

5. Host-parasite cycles when <j a is small: Compare the direction of host-parasite cycles in Figure 
11(a) to Figure 11(b): For <r Q = 0.15, both host and parasite exhibit their peaks at almost the same 
time, but when <r Q = 0.05, the host peaks at the lowest population of the parasite and the parasite peaks 
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Figure 7: The stable interior equilibrium is (0.5428,1.0841,0,0). The dynamics of the co-evolutionary model (5) when 
< Tjfj = = A'oi = r i = K 02 = 1. V 2 = 0.25 > d = 0.185, ao = 5, e = 0.9, h = 4, <r a = 0.5 with initial values 

(xi(0), 2 : 2 ( 0 ), rti(0), U2(0)) = (0.5,2,1,0.1): the host xi is red; the parasite X 2 is blue; the trait value for the host «i is green; 
and the trait value for the parasite u 2 is cyan. 


at the lowest population of the host. These simulation results suggest that the ability of the host and 
parasite to adapt, as measured by the variance of the trait differences in the parasitism efficiency a , can 
alter the community dynamics of natural systems, leading to novel dynamics including antiphase and 
cryptic cycles. Such cycles, induced by the variance of the trait difference measured in the parasitism 
efficiency a, could be another potential signature of host-parasite co-evolution and reveal that host- 
parasite co-evolution can shape, and possibly reverse, community dynamics. These results offer another 
perspective on the recent work by Cortez and Weitza (2014), which uses an eco-co-evolutionary model 
to show that predator-prey co-evolution can drive population cycles where the opposite of canonical 
Lotka-Volterra oscillations occurs; predator peaks precede prey peaks. 


5. Conclusion 

Host-parasite co-evolution, or the reciprocal evolution of host defenses and parasite counter-strategies, 
can affect a range of ecological and evolutionary processes (Thompson 2005; Gomez et al. 2014), from pop¬ 
ulation dynamics (Yoshida et al. 2003), to the maintenance of genetic variation (Clark et al 2007). Despite 
this, there are relatively few mathematical models examining the co-evolution of quantitative traits in hosts 
and parasites (Best et al 2009). Although most models assume that the parasite is obligate, cases of facul¬ 
tative and/or generalist parasites are common and ecologically significant (Spottiswoode et al 2012). 

In this paper, we have presented a co-evolutionary model of a social parasite-host system that includes (1) 
ecological dynamics that feed back into the co-evolutionary outcomes; (2) consideration of both obligate and 
facultative parasitic strategies; and (3) Holling Type II functional responses between the host and parasite. 
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Figure 8: The oscillating dynamics when rr n is small. The dynamics of the co-evolutionary model (5) when = &k 2 = Koi — 
ri = Kq 2 m 1, r 2 = 0.25 > d = 0.185, ao = 5, e = 0.9, h = 4, rr a = 0.15 with initial values (a;i(0), £2(0), rti(0), «2(0)) = 
(0.5, 2,1, 0.1): the host x i is red; the parasite x 2 is blue; the trait value for the host u i is green; and the trait value for the 
parasite U 2 is cyan. 


The analytical study on the proposed model provides insightful information on the impacts of host-parasite 
co-evolution on varied ecological and evolutionary processes. Here we summarize our main results as follows. 

Recall that d is the death rate of a parasite due to searching/hunting for all potential hosts, and r 2 > 0 
is the intrinsic growth rate of the parasite without parasitizing a given host x 2 . If d > r 2 , the parasite is 
obligate; and if d < r 2 , the parasite is facultative. In the absence of evolution, we performed local and global 
analyses to investigate the ecological outcomes when parasites range from facultative to obligate. 

When we fix other parameters and let d vary, our proposed ecological model can exhibit a wide range of 
dynamics illustrating how ecological dynamics change when the parasite makes the transition from faculta¬ 
tive to obligate. The typical dynamics are shown Figure 1-2. As an example displayed in Figure 2, we can 
see that: when d is extremely small, the host goes extinct globally; when d is small, the system can exhibit 
bi-stability between the parasite-only boundary equilibrium and the coexistence interior attractor; when d 
is in the intermediate range of values, the system is permanent and can process three interior equilibria that 
lead to two distinct coexistence attractors; when d is large, the system is permanent with only one coexis¬ 
tence attractor; however, when d is extremely large, the parasite goes extinct, resulting in global stability of 
the host-only boundary equilibrium. More specifically, our analytical results imply that: 


1. Facultative parasitism can drive the host extinct locally (Proposition 3.1) or globally (Theorem 3.3), 
while obligate parasitism can only generate global extinction of parasites (Proposition 3.1 and Theorem 
3.3). 
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Figure 9: The larger oscillating dynamics when a a is smaller. The dynamics of the co-evolutionary model (5) when 
(Jk 1 = a K 2 = K 0i — r i = K 02 = 1, V 2 = 0.25 > d = 0.185, ao = 5, e = 0.9, h = 4, cr a = 0.05 with initial values 
(a^i(0),# 2 ( 0 ), wi(0), 1 ^ 2 ( 0 )) = (0.5,2,1,0.1): the host x\ is red; the parasite #2 is blue; the trait value for host u\ is green; and 
the trait value for parasite U 2 is cyan. 


2. When the parasite is obligate, i.e., d > 7 * 2 , the host always persists. However, when the parasite is 
facultative, i.e., d < r 2 , the host can go to extinction under certain conditions, while the facultative 
parasite always persists (see Theorem 3.1). 

3. The host-facultative parasite model can have rich dynamics that process one, two, and three interior 
equilibria with consequences of two or three attractors; the host-obligate parasite model can have either 
one or three interior equilibria (Theorem 3.2). When the system has two interior equilibria (only for 
the facultative parasite) the system has a parasite-only boundary attractor and a coexistence interior 
attractor; when the system has three interior equilibria the system has two distinct coexistence interior 
attractors. 


Host-parasite co-evolution plays an essential role in both ecological and evolutionary processes. Our 
work on the co-evolutionary dynamics confirms this, and addresses the importance of trait function effects 
on evolutionary outcomes. More specifically, our main findings are: 


1. When the death rate d of the parasite depends on its trait value, under proper trait functions (Theorem 
4.1), the parasite can choose different strategies such that it can be facultative at some trait values 
while it is obligate at others. However, these strategies are not ESS (see Figure 3-4). A potential bio¬ 
logical example supporting these results may be slave-making ants in the Formica sanguinea complex. 
These species display the ability to vary behaviorally between reliance on parasitism of other ants, and 
producing their own worker offspring. They have been suggested to represent an intermediate parasitic 
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Figure 10: The oscillating dynamics converge to the stable interior equilibrium is (0.5428,1.0841,0,0) when cr a = 0.005. The 
dynamics of the co-evolutionary model (5) when ctk 1 = cjk 2 = = r\ = Kq 2 = 1, 7*2 — 0-25 > d = 0.185, ao = 5, e = 

0.9, h = 4, cr a = 0.05 with initial value (ici(0), £ 2 ( 0 ), iti(0), it2(0)) = (0.5,2,1,0.1): the host x\ is red; the parasite X 2 is blue; 
the trait value for host u\ is green; and the trait value for parasite U 2 is cyan. 


stage, between those species without predominant brood raiding, and obligate slave-making species 
(Topoff and Zimmerli 1991; Ruano et al. 2013). 


2. When the death rate d of the parasite is independent of its trait value, Corollary 4.1 implies that 
the host-only equilibrium E X1 o ulU2 can have ESS when the parasite is obligate; and the facultative 
parasite-only equilibrium E 0x2U1U2 can also have ESS. 

3. Let the death rate d of the parasite be independent of its trait value, and other trait functions follow 
normal distributions. Our results (Theorem 4.3 and its corollary 4.2; Theorem 4.4) show that: (1) 
co-evolution can save the host from extinction by destabilizing the facultative-parasite only boundary 
equilibrium and generating a new locally stable coexistence interior equilibrium (see Theorem 4.2 and 
Figure 5-10); (2) the variances of the proposed trait functions, i.e., cr^y, uk 2 and er 4 , play essential roles 
in guaranteeing the local stable interior equilibrium (ay, x%, 0,0) having ESS; (3) when the parasite 
death rate d is large, the host can have ESS that drive the parasite extinct globally; and (4) Large 
variances of the parasite carrying capacity <jk 2 and of the parasitism efficiency a a are required to make 
sure of (a^, ^ 2 ,0,0) being ESS. In addition, we have following interesting results: 

(a) It is possible to have locally stable coexistence interior equilibrium (ay, x%, it*, 11 %) with u*, > 0 

(see Proposition (4.1) and Figure 5-6) but these strategies may not be ESS. This can occur when 
the ecological system has two interior equilibria with the facultative parasite-only equilibrium 
being locally stable; and co-evolution destabilizes this facultative-parasite-only boundary equilib- 
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rium and generates a new locally stable coexistence interior equilibrium . 


(b) The variance of the trait difference in parasitism efficiency has huge impacts on the population 
dynamics. More specifically, small values of the variance cr a can destabilize the system, thus 
generate evolutionary arms-race dynamics with different host-parasite fluctuating patterns (see 
Figure 5 - 11). 


Our theoretical work on the ecological and co-evolutionary dynamics of the host-parasite system show 
interesting matches with recent ecological considerations of social parasitism (Foitzik et al. 2001 & 2003; 
reviewed by Kruger 2007; Davies 2011; Kilner and Langmore 2011) in the following ways: 

1. The study in Section 4.1 of the host-only equilibrium suggests that the host-only equilibrium E xi0uiU2 
of the co-evolutionary model (5) can have two ESS (u^u^,) = (0, c) and (u^u^) = (c, 0) when the 
parasite is obligate (i.e., d > V 2 ). This can be classified as one of the co-evolutionary outcomes when 
the host successfully resists invasion by the parasite, and when resistance is due to effective front-line 
defenses (Foitzik et al. 2001 & 2003; reviewed by Kruger 2007; Davies 2011; Kilner and Langmore 
2011 ). 

2. The results in Theorem 4.2, Theorem 4.3 and its corollary 4.2 imply that large variances of the parasite 
carrying capacity uk 2 and of parasitism efficiency <r 0 are required to make sure of the strategy ujjl) = 
(0, 0) of (x\, ^ 2 ,0,0) being an ESS but it may not be a global ESS, because of the possibility that the 
system has other ESS for coexistence. This can be considered as equivalent to a tolerance of the parasite 
co-evolutionary outcome where hosts to some degree concede to the parasite, and accept it into their 
nests, but also make adjustments to their life history (or other traits) to minimize the negative effects 
of parasitism on their fitness (Foitzik et al. 2001 & 2003; reviewed by Kruger 2007; Davies 2011; Kilner 
and Langmore 2011). This co-evolutionary outcome is more likely to be the case that parasitic control 
of the co-evolutionary trajectory. 

3. Theorem 4.4: A large death rate due to the parasite overhunting/attacking the host can lead to the 
extinction of the parasite, and lead to a successful resistance by hosts evolutionary outcome (Foitzik 
et al. 2001 & 2003; reviewed by Kruger 2007; Davies 2011; Kilner and Langmore 2011). 

4. Theorem 4.4: when the ratio of the variance of the parasite to the host is larger than 1, i.e., > 1, 

small values of ag can lead to the interior equilibrium (ccl,a^, 0,0) being the only equilibrium that has 
ESS in the co-evolutionary model (13). This type co-evolutionary outcome is referred to as acceptance 
of the parasite. It can become an adaptive strategy for hosts when the costs of rearing a parasite are, 
on average, lower than the costs of recognition and rejection (Foitzik et al. 2001 & 2003; reviewed by 
Kruger 2007; Davies 2011; Kilner and Langmore 2011). 


6. Proofs 

Proof of Lemma 3.1 

Proof. Both sets X\ and X 2 are invariant for the ecological model (6) since we have 

dxi I _ dx 2 I _ 

~df\xi=o ~ ~fff\x 2 =0 ~ 

The continuity of the expressions of both and implies that Model (6) is positively invariant in X. 
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Let Ki ' 1 = sup UieR {A'i(wi)} and K ^ 1 = sup U2gR {A 2 (it 2 )}. From (6) and its positively invariant property 
in A, we have the following two inequalities 


^ = *r rrU-ft)-^ < n* l( l - &) < n» l( l - 


K i - 


(fa? 

dt 


= #2 


TS _-d + r 2 (l-^) <f+r 2 , 2 (l-^) = (r2 + fh 1- 


K 2 M (-2 + f) 


< [r 2 + f )a; 2 1 - 


K 2 U ^2 + -f) 


which implies that 


limsupa;i(i) < and lim sup x 2 (t) < 

t—too t—>oo 


K^{t 2 + f) 
r 2 


Therefore, the statement of Lemma 3.1 holds. 


□ 


Proof of Proposition 3.1 

Proof. The Jacobian matrix of the ecological model (6) evaluated at its equilibrium can be repre¬ 

sented as follows: 


J(xl,x*) 


Therefore, we can have the following three cases: 

1. The eigenvalues of the extinction equilibrium Aoo are Ai = rq and \ 2 = r 2 — d. Thus, Eqq is a source 
if r 2 > d while it is a saddle if r 2 < d. 

2. The eigenvalues of E w are \\ = —rq and A 2 = r 2 — d+ - [ ™ 0 ^ 1 Kl ■ Thus, E w is a sink if r 2 + < d 

while it is a saddle if r 2 + 1 + a 0 %} <1 > d. 

3. The equilibrium Aoi exists if r 2 > d. In this case, Aoo is a source and Aio is a saddle. The eigenvalues 

of A 0 i are Ai = rq + ad Jf 2 — aK 2 and \ 2 = d — r 2 <0. Therefore, Aoi is a sink if rq + adI ^ 2 < aK 2 
(i- e -: ^7 < 1 - while it is a saddle if rq + > al < 2 (i.e., > 1 - ^). 

Thus, the statement of Proposition 3.1 holds. 

□ 


/ 1 2'1 \ (LX n . 

m 1 - K-J ~ T+ht? + 


* ( nx£ 

c i v — 


(1 -\-ahx^) 2 ^ 


1 -\-hax^ 


(1 -\~ahxi ) 




K 2 1 -\-hax^ 


)-d-^ 


( 21 ) 


Proof of Theorem 3.1 


Proof. According to Lemma 3.1 and Proposition 3.1, we can conclude that (i) Model (6) is positively invariant 
in both Ai and A 2 ; (ii) the omega limit set of Ai is Ado; and (iii) the omega limit set of A 2 is Eqo if r 2 < d 
while the omega limit set of A 2 is A 01 if r 2 > d. 

The results (Theorem 2.5 and its corollary) in Hutson (1984) guarantee that the persistence of the host 
X\ is determined by the sign of 


daq 

X\dt 


x 1 


dx 1 
x\dt 


if V 2 < d 

Eqo 


dx 1 
xi dt 


Eqi 


if r 2 > d 
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And the persistence of the parasite a ; 2 is determined by the sign of 


dX 2 


X2 dt 


dx 2 


x. 


X2 dt 


E 1{ 


Therefore, we can conclude the following: 

1. The prey X\ is persistent in X if the following inequality holds 


dx 1 


X\dt 


dx 1 


a' 2 


x\dt 


E 10 


dx 1 
x\dt 


dx 1 
x\dt 


En, 


Eqi 


= ri > 0 if r 2 < d 


(22) 


= n - o(l - ~ )K‘2 > 0 if r 2 > d 


This implies that X\ is persistent if r 2 < d or > 1 — > 0. 

2. The predator X 2 is persistent in X if the following inequality holds 


dx 2 


X2 dt 


dx 2 


Xx 


X2dt 


, eaK[ 

= r 2 — d + --—— > 0 

E 10 1 + riaKi 


(23) 


This implies that X 2 is persistent if r 2 + > d. 


3. The discussions of two items above imply that Model ( 6 ) is permanent (i.e., both prey X\ and predator 
X 2 are persistent in X) if either 


r 1 d 

> 1 - >0 


or 


T 2 


aK 2 r 2 

eaK 1 


1 + haK\ 


> d> r 2 . 


□ 


Proof of Theorem 3.2 

Proof. The interior equilibrium ( x\,X 2 ) of Model ( 6 ) satisfies the following two equations: 


ri(l 


XX \ 
Kx> 


ax 2 

l-\-haxi 


0X2 = f(x 1 ) 


rx(Kx- Xl )(l+ah Xl ) subject to Q < ^ 


eax 1 
l-\-hax\ 


d-r 2 ( 1 - 
subject to 


k 2 


j U- 2 








< *1 and d < r 2 + 


(24) 


which implies that there is no interior equilibrium {x\,x^) if d > r 2 + ■ The equations (24) imply 

that x\ is a root of fg{x 1 ) = g{x 1 ) — f(x\): 


fg(xi) = 


C3X1 + C 2 x\ + C\X\ + cO 
ar 2 Ki(l + ahx\) 


F {xi) 

ar 2 A'i(l + ahxi) 


where 


C 3 = rir 2 (a/i ) 2 c 2 = a/irir 2 (2 — ahK 1 ) 

ci = a 2 hl<il <2 (r 2 + | - d) + rqr^l - 2ahK 1 ) c 0 = ar 2 A'i A ' 2 ^1 - ^0 

We classify interior equilibria for Model ( 6 ) in the following cases: 


n 

a 
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1. If Co > 0 <t=> K 2 ^1 — ^ > 0 <t=> < 1 — ^ (i.e., sufficient conditions for the locally asymptotical 

stability of Eqi according to Proposition 3.1) , then we can conclude that F(x i) has either no interior 
equilibrium or two interior equilibria since C3 > 0. This implies that Model (6) may have either no 
interior equilibrium or two interior equilibria when Co > 0 (i.e., Eq 1 is locally asymptotically stable). 

2. If Co < 0 <t=> > 1 — ^ (i.e., E 01 is a saddle according to Proposition 3.1), by using similar arguments 

for the case of Co > 0, then we can conclude that Model (6) may have either one interior equilibrium 
or three interior equilibria when cq < 0. 

3. Note that 

F'{x 1 ) = 3c 3 xl + 2c 2 Xi + ci = 3c 3 (aq + + Ci - ^ 

= 3c 3 (Si - ^l ) 2 + a 2 r 2 K x K 2 (h + ± - ™(o/.* 1+ i) 2 

Thus, we have the following two scenarios: 

(a) If a 2 r 2 KiK 2 (h + ^ — -^j < rir2 ( Q ^i+ 1 ) (i.e., Cl < ), then F'(x 1 ) = 0 has two real roots 

1 — C2 — a/ C 2 — 3ciC3 2 —C2 + a/ c 2 ~ 3ciC 3 

x c =- \ - <x c = --■ 

3C3 3C3 

In this case, we can conclude that F(x 1 ) has three positive roots if 

c 0 < 0, x\ > 0, F{x\) > 0 and F{x 2 ) < 0. 


The equation F{x 1 ) has one positive roots if 

co < 0, and F(x 2 ) < 0. 

The equation F(x 1 ) has two positive roots if 

c 0 > 0, x 2 > 0, and F(x 2 ) < 0. 

(b) If a 2 r 2 KiK 2 (ji + £ — > r ^ r ^ ah ^ +1 ) (i. e ., Cl > S 2 .), then F'{x 1 ) > 0 which indicates that 

F[x 1 ) is an increasing function of X\ . Therefore, F(x 1 ) has one positive root if co < 0 and has no 
positive roots if co >0. 


4. Now we discuss that the cases when Model (6) has either no interior equilibrium or one interior 
equilibrium in terms of the function f(x 1 ) and g{x 1 ). First, we know that (1) The function f(xi) is 
a degree two polynomial with two roots — ^ and K 3 . (2) f(x\) is a decreasing function whenever 
X\ > Kl ~^/ ah w hHe g{x\) is an increasing function with a unique positive root x* = a r e J^l r2 ) ] when 
r 2 < d < r 2 + f. Therefore, the functions f(x±) and g(x 1 ) have no intercept for X\ > 0 if 


d— r 2 

a[e - h(d - r 2 )\ 


> K t <=> d > r 2 + 


eaKi 
1 + ahKi 


This also implies that when r 2 < d < r 2 + we have the following sufficient condition for f(x 1 ) and 
g(x 1 ) having a unique positive intercept: 


d — r 2 


< K 1 r 2 < d < r 2 + 


eaKi 


a[e — h(d—r 2 )\ ' ^ ' 1 + ahKi 

If r 2 > d, then the functions f(x 1 ) and g(xi) have no intercept for xi > 0 in the following two cases: 
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(a) < o and g( 0) > /(0) ^ K 2 (l - £) > 

(b) > 0 and ff ( 0 ) > /(^*) O tf 2 (l - 4 ) > ri( ^°'g l)2 - 

Now assume that K\ < -K-, then f(x i) is a decreasing function while g(x\) is an increasing function 
for X\ > 0. Thus the functions f(x i) and g(x i) have an unique interior intercept for X\ > 0 in the 
following two cases: 

(a) r 2 > d and g( 0) < /(0) <;=> 0 < A" 2 (l - 

(b) r 2 <d<r 2 + f and 0 < x* = a [ e ^h(d-r 2 ) ] < where = °- 

The discussions above lead to the following summary regarding the interior equilibria of Model (6) which is 
also listed in Table 2: 


1. Model ( 6 ) has no interior equilibrium if one of the following conditions is satisfied: 
( a ) rf>r 2 + or 


(b) ^< 1 - 


I d _e_ 

' r 2 r 2 i 


r 1 r 2 (ahK 1 + l) 2 . 

3 ’ ui 


and a 2 r 2 K\K 2 (h 

(c) /f 1 <-Land^<l-A ;0 r 

(d) ^>^and r f+^ 2 <l-^. 

2. Model ( 6 ) has an unique interior equilibrium if one of the following conditions is satisfied: 
(a) r 2 <d<r 2 + or 


(b) 


aK 2 


> l- 


-d- and (h + — — 

r2 \ r 2 r 2 J 


> 


ri (ahKi~\~l) 2 # 


3a 2 K 1 K 2 '■ 0r 

(<=) ^ > 1 ~ (fe + 77 ~ 4r) < and F(xl) < 0. 


3. Model ( 6 ) has two interior equilibria if < 1 — and F{x 2 ) < 0 with 

x\ > 0. 

4. Model ( 6 ) has three interior equilibria if ^ > 1 - (h + £ - < ri{ 3 , ^(^c) > and 

F(x 2 ) < 0 with x\ > 0. 


where a* = ^ - y^cT 


-3_ < a;2 _ C2_ 

3c 3 ^ "‘'c 3 c 3 


ci - 


3 c 3 ' 


Let be an interior equilibrium of Model ( 6 ), then its stability is determined by the eigenvalues of 

the following Jacobian matrix 


J, 




r * ( r l I ea:r 2 A 


X 1 V + (1+a/dq)^ J 

l+/iax* 

eax 2 

r 2 x^ 

(1+a/irc^) 2 

K 2 


(25) 


whose eigenvalues A,, i = 1, 2 satisfy the following equalities: 

riX* r 2 x 2 ha 2 x\x 2 


Ai + A 2 — — ■ 


A', 


and 


AiA 2 = x\x* 2 


ri r 2 


K 2 (1 + ahx {) 2 
r 2 a 2 hx I 


< - 


rqx | r 2 x 2 ax 2 


K , 


Ab 


KiK 2 (1 + a/ix *) 3 (1 + a/ix *) 2 


> * * 
_ x x x 2 


ri r 2 


KiK 2 (1 + ahKiY 


- a 2 hK 2 (r 2 + -) 
h 
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This implies that the interior equilibrium (a;!, 22 ) is locally asymptotically stable if 


ri a rir 2 

ir 2 > 2 “ d Kjr, 


ea 


(1 + ahK i) ; 


> a 2 h.K 2 {r 2 + j). 


Therefore, the statement of Theorem 3.2 holds. 


□ 


Proof of Theorem 3.3 

Proof. The proof of the global stability of E m is similar to the proof of the global stability of E w , thus we 
only focus on the case of £01 • 


According to Proposition 3.1, we know that Model ( 6 ) has only three boundary equilibria Eqq, E\q and 
E 01 where £00 is always a saddle; E w is locally asymptotically stable when the inequality d > r 2 + i^hkx 
hold while it is unstable if d > r 2 + ; and £01 is locally asymptotically stable when the inequality 

< 1 — ^ hold while it is unstable if > 1 — Proposition 3.1 also implies that if £10 is locally 
asymptotically stable then £qi is unstable while £qi is locally asymptotically stable then Eiq is unstable. 


Now assume that < 1 — -A and one of the following conditions is satisfied 


L. a 2 r 2 M 2 (h+£-^) 


rir 2 (Q/iiCi + l ) 2 . 

^ 3 • U1 


2 - K i < ot 


3 - ^ and -f+t k 2 )2 < 1 ~£- 


Then according to the proof of Theorem 3.2, Model ( 6 ) has only three boundary equilibria £ 00 , Eio and Eqi 
where both Eqq, Eiq are unstable and Eqi is locally asymptotically stable. 

According to Lemma 3.1, Model ( 6 ) has a compact global attractor. Thus, from an application of the 
Poincare-Bendixson theorem (Guckenheimer and Holmes 1983 [39]) we conclude that the trajectory starting 
at any initial condition living in the interior of R+ converges to one of the three boundary equilibria Eqq, E\q 
since Model ( 6 ) has no interior equilibrium under the condition. Since Eqi is the only locally asymptotically 
stable boundary equilibrium, therefore, every trajectory of Model ( 6 ) converges to Eqi for any initial condition 
taken in the interior of R+. This implies that Eqi is global stable. Therefore, our statement holds. 

□ 


Proof of Theorem f.l 
Proof. Let 


a v 


da(u, v) 
dv 


a u 


da(u,v) d 2 a(u,v ) d 2 a(u,v ) A d 2 a(u,v) d 2 ci(u,v) 

- , ci vv = - 5 O'uu = 9 ' and. CL UV — “ “ = “ ~ 

au av z au z avau auav 


Now we are exploring the local stability of E Xl0uv = (AA(u), 0, u, v) and E 0x2 uv — (o,K 2 {v) (l - ,u,v 
where u, v are the trait values that satisfy the following two equations: 


r\x\K x ( u ) 
K'i(u) 


( l-\-ha(u,v)xi ) 2 


X 2 d u ( u,v ) 
(l-\-ha(u,v)xi )' 2 


d'{v) 


exi a v (u,v) , T2X2K 2 {v) 
( l+ha{u,v)x \)‘ 2 ' K£( v ) 
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Therefore, at the equilibrium E Xi q uv = (A'i(u), 0, u, u), we have 


K[(u) = 0 and d'(v) = 


eKi{u)a v {u, v) 

(1 + ha(u, v)Ki(u )) 1 2 


with the following Jacobian matrix 


Jt 


(K i (u), 0 ,u,v) 


~r i 
0 
0 


K i ( u)a(u,v ) 


l-\-ha{u,v)K\ ( u ) 

e aQu,-q)Ai(u) _ ,, n 

l+ha(u,w)ii'i(u) “W + f2 
< 7 ? a u (u,v) 
(l+ha(u,v)Ki ( u )) 2 

a%ea v (u,v) (1 — ha(u,v)K\ (u)) 

(l+/ia(u,u)-Ki(u )) 3 A|(t>) 


0 

0 

a?riA'"(i 

A'i(u) 

&21 


0 

0 

0 

b 2 2 


(26) 


with 


and 


&21 = 


&22 = 


< 7 2 eKi{u) [a uv (u, u)(l + haK\{u )) — 2Ki{u)a u {u, v)a v (u , u)] 
(1 + ha(u, v)Ki(u )) 3 

a 2 eKi(u) [a vv (u, u)(l + haK\{u)) — 2Ki(u)a%(u,v)] 


(1 + ha(u, v)Ki(u )) 3 


- °ld"(v). 


The eigenvalues of (26) are 
ea(u , v)Ki(u) 


rfriK"(u) cr%eKi(u) [a vv (u,v){ 1 + hciK 1 (u )) - 2K 1 (u)a%(u, u)] 


~ ri ’ 1 , 7 , ' \rf / . -d(^)+r 2 , — , . , -- x,, , - v 2 d W- 

1 + ha(u, v)Ki(u) Ai(«) (1 + lia{u, v)hi(u)y 

This implies that if the equilibrium E Xl q uv = (ATi(u), 0, it, u) exists, then it is locally asymptotically stable if 

ea(u,v)K 1 (u) eii'i(u) [a vv (u,v)(l + haK^u)) - 2i^ 1 (u)a^(u,u)] , 

1 + ha(u,v)K 1 (u) w (1 + ha(u,v)K 1 (u )) 3 w 

Similarly, if the equalities o u (u, v) = 0, d(v) < r 2 and d'(v) = r2 ^ 2 2 ff hold, then the equilibrium Eq X2UV = 

^0, K 2 (u) ^1 — , u, vj exists with the following Jacobian matrix 


Jt 


(0,K2(v),U,v) 


where 


n - a{u,v)K 2 (v) (l - 0 

ea(u,u)Ar 2 (v) (l - ^) d(v) - r 2 


CTjriifJ(u) 

Ajf(«) 

CT 2 ea.^, (zi, u) 


0 -ofA^v) (l - a„„(u,u) 6 


K^v) 


d{v) 


0 

0 

12 

&22 


(27) 


612 = -<t(K 2 (v) ( 1 -— ) a uv (u, v) 


and 


&22 = 


^2 (1 - [A' 2 (u)A^(u) - 2{K’ 2 {v)) 2 } 


K 2 {vf 


- a 2 d"(v). 


The eigenvalues of (27) are 


r\ — a(u, v)K 2 {v) 



~d{v) +r 2 , 


-alK 2 {v) 



(^b ^) ; 
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and 


° 2 r 2 (l - tt) [ k 2 (v)K%(v) - 2 {K' 2 (v)Y 


- °ld'\v). 


K 2 {vf 

This implies that the facultative-social-parasite-only equilibrium Eq X2UV = (0, K 2 (v), u, v) exists and it is 
locally asymptotically stable if the following inequalities hold 


r i 


1 - 


d{v 


< a(u,v)K 2 (v), r 2 < d(v), a uu (u,v) > 0, and 


r 2 (l - ^) [K 2 {v)K'f{v) - 2 (K! 2 (v)Y- 
K 2 (v) 2 


< d"(v). 

□ 


Proof of Theorem f.2 

Proof. First, we can check that the boundary equilibrium of (13) can occur only if u\ = u 2 = 0. According 
to Proposition 3.1, the boundary equilibria are: 


Eqooo — (0,0,0,0), A’xiooo — (Aoi,0,0,0), 


and EqxzOO = 



provided that r 2 > d. 


When u\ = u 2 = 0, the Jacobian matrix evaluated at an equilibrium {x\, 0, 0) can be represented as 


follows: 


Oil 

Ol 2 

0 

0 



J(x^ ,x% , 0 , 0 ) 

a 2 i 

0 

a 22 

0 

0 

bn 

0 

bi 2 

(28) 



0 

0 

b 2 \ 

b 22 



whose eigenvalues are determined by the eigenvalues of matirices A and B 


eapx ^ 
l-\-aphx^ 

(29) 

and 

CLQX 2 \ CT-^CLpX^ 

{l+aphx^) 2 J cr 2 (1-j-aohx ^) 2 

_ 2 ( r 2 x% ,_ cipexl 

fF 2 y a K 2 k 02 a'i{.l+a Q hxiy 2 

Therefore, we can conclude the following cases: 


B = 


hn bi 2 
b 2 i b 22 


r ii, 
r Ki K oi 


+ 


<t 2 eapx j 
a 2 (1+aohx 



A = 


a n ai 2 
«21 022 


* ( n i ea 0 x ? \ , ( /l _ x l \_ apx? 1 

1 ^ K 01 (l+a 0 hxl)' 2 J ( ^ iCoi 1 +a 0 hxl J 


l-\-hapx £ 


eapx 2 

(l+a a hxl ) 2 * 




1. -Eoooo : In this case, the eigenvalues of A are r\ and r 2 — d while the eigenvalues of B are zeros. Thus, 
Aoooo is always a saddle. 

2. E Xl ooo = (/ioi, 0,0, 0): In this case, the eigenvalues of A are — n and r 2 — d + while the 

2 2 TV- 

eigenvalues of B are — and — (i+aohK 01 y 2 • Thus, A’xiooo is locally asymptotically stable if r 2 + 

1 +ahKm < ^ (i- e -> ^ 10 l° ca Hy asymptotically stable for the ecological model ( 6 )) and i? xl ooo is a 

saddled r 2 + I ^ r >d. 

3. E 0 x 2 q 0 = (0, K 02 ^1 — , 0,0) exists if r 2 > d: In this case, the eigenvalues of A are r 1 + a 0 ^ 02 — a 0 K 02 

and d — r 2 while the eigenvalues of B are — °~i ho 20 p(d-r 2 ) anc [ a 2 (d-i 2 ) . t^us, E 0x 00 is always a saddle. 

°'a r 2 a K 2 
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4. When = uZ, = 0, the existence conditions regarding interior equilibria (x*,x 2 ,0, 0) of the evolu¬ 
tionary model (13) are the same as interior equilibria (x*, x 2 ) of the ecological model (6) by letting 
Ki = A'oi. K 2 = K 02 and a = a 0 . Thus, according to Theorem 3.2, Model (13) can have one, two, or 
three interior equilibria (x*, x 2 , 0, 0) depending on the values of ri,d,h,e,ao and K<h where i = 1,2. 
Sufficient condition of an interior equilibrium (x^, x|, 0,0) being locally asymptotically stable is that 
both matrices A and B have negative eigenvalues. 

According to Theorem 3.2, the eigenvalues of A are negative if inequalities (7) hold, i.e., 


^2 Uo . n r 2 

K02 2 and K 01 Ko 2 


(1 + a 0 hK 01 ) 3 


> Oq hK 02 (r 2 + -). 


Let Ai and A 2 be eigenvalues of B. Then we have as follows: 


Ai + A 2 = — 


rifj X 1 

K 01 °K 1 


< - 


K 02 v'k 2 


ao(<r^x^ — ea 2 x 2 ) 
( l+a 0 hxi )' 2 




- ao 


2 ^ T 2 (J 2 ^2 

~ K 0 2 *^ 


aoecr 2 x 2 
( 1 +a 0 hxl ) 2 


and 


A 1 A 2 — 


'l u 2 


> 


' l u 2 


rir 2 x 1 x 2 


K 01 K 0 


r 1 r 2 x 1 x 2 

Kq\Kq 2 


riV 2 X 1 x 2 


r\r 2 x^ x 2 


aor 1 e(xl ) 2 


a K 1 Koio-2(l-\-ahxl) 2 
r 2( x 2 ) 2 


aor 2 (x 2 ) 2 


a K 2 K 02 < j2 {l-\-ahx *[) 2 


rje(xp 2 

K 01 




0 - 2(1 -fa/ix *) 2 


_ aorie(x *) 2 _ CLQr 2 {x 2 y 

KoiK 02 1 o-^ i K 01 ( 7 %(1+ahxl ) 2 a‘j <2 Kq 2 (t^ l 


KqiKq 2 


apr 2 (x 2 ) 2 , a 0 rie(x*) 2 

ct‘ 2 K2 Kq 2 g 2 K 01 a% ( 1+ahxl ) 2 


T 2 X 2 I nx-y 
K 02 [ K 01 


a ox 2 
0-2 0-2 

k 2 a 


4- 


qprie(x*) 2 


a K 1 K 0l&a (1+ahxl ) 2 


Therefore, we can conclude that if K J 1 1 a i > and one of the following two inequalities hold 


*1 


n_Xj 

K 01 


> 


riecr 2 K 2 K 02 (, 

r 2 <?K 1 Koi 


7^ C 
K 2 


^ r l (7 K 2 C7 a . X* 

<t=> -; tA - > -Jr Or 


aoKoi 

_ 2 2-k Gk 2 r\eKp 2 ^ (x 2 ) 2 

W) 2 ffjfj r 2 K 01 > (*f) 2 


(31) 


then we have Ai + A 2 < 0 and A 1 A 2 > 0, which implies that B has two negative eigenvalues. Thus, we 
can conclude that if the interior equilibrium (x^,X 2 , 0 , 0) exists, then it is locally asymptotical stable 
if one of the two inequalities in (31) holds and the following inequalities hold: 


f 2 r 2 n -i rir 2 

1 ~P —> 2 ir ( > °° and ts + 

i K(yl Aqi J Aqi A 02 


GCL n 9 1 T/r , G V 

(i+^,,p > “»'* Ao2( ’' 2 + 


Therefore, the statement of Theorem 4.2 holds. 


□ 


Proof of Theorem 4.3 

Proof. According to the definition of ESS, A^ooo is an ESS if it is an ESE and it satisfies the ESS maximum 
principle (2). According to Theorem 4.2, E X1 000 is an ESE if the inequality r 2 + d holds. Notice 

that 


max„ l6 R{Gi(ui, 0,0, Aqi, 0)} = max t , l€R jn 


K 0 


= ri (l - 


0} 


K oi 

K\{v\) 


Ki{vi) t 

= G\ (0,0,0, Af 01 ,0) = 0 


vi =0 
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and 


max,, 2 g R {G 2 (u 2 ,0,0, A' O i,0)} = max U 2 g R j ~d + r 2} 


e a( 0 ,v 2 )K 01 
1 +h a( 0 ,V 2 )Koi 


- d + r 2 


V 2 —0 


= G 2 (0,0,0, K 0 i, 0) < 0 


Therefore, E X100 o is an ESS whenever the inequality r 2 + < d holds. 


Similarly, the interior equilibrium (x*, x 2 , 0,0) is an ESS if it is an ESE and it satisfies the ESS maximum 
principle (2). According to Theorem 4.2, (x*,x 2 ,0, 0) is an ESE if the inequalities (18) hold. Notice that 
Kiiv ), a( 0 , v), a(v , 0 ), i = 1 , 2 are decreasing functions with respect to |u|, and the following equations 


G 2 (v 2 , 0 , 0 , a:*, *2) 

dCr2 (^2 , 0 , 0 , 21 * , 2?2 ) 

dv 2 


Therefore, we can conclude that 


max U 2 eR {G 2 (u 2 ,0,0, x\,xl)} 


On the other hand, we have 
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- d + r 2 i 1 - TW)} 
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Since the inequalities (18) hold, then we have 


iQi.Q) 
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max I 1^2. l r 1 eK 02 a K 2 ^ \ 
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‘ \ a Ki V r 2 K 01 ’ a 0 K 01 J 


-, then we have 


aoi(«,, 0,0.*;, < 0 it m > 0 . aG 1 (»„ o.o. a; ;.x;) > „ if n < 0 






Therefore, we can conclude that 
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if the following inequality holds 
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Proof of Theorem 4-4 

Proof. According to Theorem 4.3, we know that the boundary equilibrium E X1 poo is an ESS of the co- 
evolutionary host-parasite model (13) whenever the inequality r 2 + < d holds. Under this condition, 

Model (13) has only two boundary equilibria Aoooo and E x1 qqq where Aoooo is always a saddle. 

According to Lemma 3.1, we know that limsup^^ Xi(t) < Kpi . Since a(iti,it 2 ) < a o and the function 
l+hatuiU^X! increases with respect to aaq, thus for time large enough, we have 


dx 2 
X 2 dt 


ea(ui, u 2 )xi 
1 + ha(u\,U2)x\ 


d + r 2 { 1 


x 2 

K 2 {u 2 ) 


)< 


eapKoi 
1 + ciphKpi 


d + r 2 < 0 . 


This implies that Model (13) has no interior equilibrium. Therefore, when the inequality r 2 + i^ahKoi < ^ 
holds, Model (13) has only a unique ESE E Xl00 o which is also an ESS. 

If Kp 2 ap < an< ^ ^oia 0 < ^ holds, then for any given trait u = {ui,u 2 ), we have K 2 a < 1 f 1 j_ ■ 

t '2 r 2 

According to Proposition 3.1 and Theorem 3.2, the ecological model ( 6 ) has a unique interior equilibrium 
(a;*, £ 2 ) and three boundary equilibria Aoo, E 10 , Ep\. Thus the evolutionary model (13) has at least one 
interior equilibrium (£*, xt,, 0 , 0 ) and three boundary equilibria -Eoooo> E xl ooo, Eo X2 oo which are all saddles. 
According to Theorem 4.3, the interior equilibrium (a;*, x%, 0, 0) is an ESS if (20) holds. 


Now we should show that the evolutionary model (13) has no interior equilibrium {x\,x 2 iu\,u 2 ) with 
u\u^ > 0 when (20) holds. Assume that this is not true, then let (a;^, x 2 , u*, u%) be the interior equilibrium, 
then it satisfies the equations(16) and (17) which gives the following equality: 


r-ixlulo- 2 K 2 K 2 {ul) _ x 2 
r 2 X 2 U 2 <x 2 Ki Ki{u\) ex\ 


I ena 2 K2 ulK 2 (ul) ^ = 
r 2 aj (i u^K 1 (ul) Xl a Kl \ 


K > 2 ( <"i ) 2 

eriu\Kp 2 e 2<T K 1 * 

- Xi . 


r 2 « 2 -Koi 


(32) 


In addition, the equations(16) and (17) requires u\,u 2 and u\ — U 2 having the same sign. This implies that 
l u il > 1 ^ 2 1. Thus, according to (32), if ctk 2 > ok x , then we have 


(T k 2 


A 

~\ 


(“|> 2 j K > 2 

er\u\Kp 2 e 2 ° K ' * a K ■ 

- Th? - x i > 

r 2 u 2 K 01 &K 


n 


K > 2 ( K > 2 

enK 02 e * (Jk 2 

- t? - *1 > - 

T 2 Eq 1 OKi 


eriAo 2 
^A'oi ‘ 


This is a contradiction to the assumption that > yv i n (20). Thus, Model (13) has no interior 

equilibrium (a;*, £’ 2 , uj, u^) with > 0 when (20) holds. This indicates that the interior equilibrium 

(£*,£ 2 , 0 , 0 ) is the unique ESS if (20) holds. 


□ 


Proof of Proposition 4-1 
Proof. If u\u *2 ^ 0, 

I eria 2 ^2 (112) 

then from equations (16) and (17), we must have £^£2 > 0. Let g{u\ 1 U 2 ) = y ’ then 

according to the equation (32), i.e. 


ri£iUicr|- 2 A'2(u2) _ £ 2 
r 2 X 2 U 2 aj Ci Ki(ui) ex\ 


o £2 = 


ena 2 K uiA 2 (u 2 ) , , 

- 2 - ——-£1 = g(u 1 ,u 2 )x 1 

£ 2 CTk u 2 Ai(u 1) 


which gives the following equality when it combines with the equation (14): 


£1 



£1 \ 


a(ui,u 2 ) 


er 1 u 1 cr] t . K 2 (u 2 ) 

2 jy / \ X 1 

r2u 2 <7| Ci K\{U\) 


1 + ha(u\,U2)x\ 


a{ui,U2)g{u 1 ,U2)xi 
1 + ha(ui,u 2 )xi 
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Thus we have a unique solution x\ that is a function of traits Ui,i = 1,2: 


= fi(ui,u 2 ) = ^ ~ 


l +aK 1 q | y/ (ri+aKig) 2 +riahKi(riahKi+2ri—2Kig) 


2 r\ah 


2 r\ah 


(33) 


where g = g{u\,u 2 ) = ■ Similarly, Equation (32) combines with (15) gives the following 

equality: 


r 2 1 - 


xig(u 1 ,u 2 ) 
K 2 (U 2 ) 


ea(ui,u 2 )xi 
1 + ha(ui , u 2 )xi 


which also gives a unique solution x\ that is a function of traits Ui,i = 1,2: 


x* = / 2 (wi,m 2 ) = ^ + 


2g 2 r- 2 ,ahg 

Therefore, we are able to solve for u\ and by letting 

r 1 u 1 


aeK 2 _ 1 _ _|_ y/(r2g—aeK2) 2 +r2ahK2(r 2 ahK 2 +2r2g+2eaK 2 ) 

2 ah ' 2 r±ah 


(34) 


fl(u*,U 2 ) = / 2 «,W2) and 


__ ( m * - U2)g{ul,U2)a(ul,u%) 

a K 1 Ki(ul) al (1 + ha{u* 1 ,u* 2 )fi{u* ll ul)y 


Based on the discussion above, we can conclude that the necessary condition for (ul, itjji, x h ^ 2 ) being an 
interior equilibrium of Model (13) is that the following equalities hold: 


*l = /l(«l,« 2 ) = 


5K)W 2 ) 


= / 2 (u*,U 2 ) and 


nui 


_ = (ut - ul)g{ul,U2)a(ul,u* 2 ) 

a Ki K i( u i) (Ja (! + ha{ul,ul)fi(ul,ul)f 
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(a) The oscillations when a a = 0.15 



t 


(b) The larger and faster oscillations when cr a = 0.05 


Figure 11: The population dynamics of the host and parasite for the co-evolutionary model (5) when ctki = &k 2 — -Koi — 
ri = Kq 2 = 1, ^2 = 0.25 > d = 0.185, ao = 5, e = 0.9, h = 4 with initial values (a^i(0), #2(0), ui(0), 7 x 2 ( 0 )) = (0.5, 2,1, 0.1), 
while the value of <j a is 0.15 (Figure 11(a)) and 0.05 (Figure 11(b)): the host x\ is red; and the parasite X 2 is blue 
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